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Abstract: In unitarity cut method, compact input of on-shell tree level amplitudes is crucial to simplify 
calculations. Although BCFW on-shell recursion relation gives very compact tree level amplitudes, they 
usually contain spurious poles. In this paper, we present a method to deal with this issue and provide 
explicit simple algebraic functions for various coefficients in the presence of spurious poles. As an applica- 
tion, we present analytic result (not just rational term) for one-loop five gluon A{ — h + + +) with scalar 
propagator for the first time. 
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1. Introduction 

The computation of scattering amplitudes at the next-to-leading order (NLO) and even beyond is not only 
a pure theoretical interest, but also a definite requirement of experiments, such as the current Tevatron 
collider and the forthcoming Large Hadron Collider (LHC) experiments. In last few years, there have been 
spectacular progresses for this subject, which are extensively reviewed in the report 

The art of the NLO calculation is based on the knowledge that, any one-loop amplitudes can be 
expanded into a set of master integrals as [§, ||, [|] 

A = c M (6)/ 5 D ' (i) [l] + c 4 , 4 (6)/ 4 D ' (i) [l] + C3,(6)/ 3 D ' (i) [l] + c 2 ,(e)/ 2 D ' (i) [l] + Cl>i (e)/f ,W [1], (1.1) 

where the master integral (i.e. scalar integral) is defined as 1 

pDfil = -iUit) D l 2 [ JlZ - fl 2) 

J i^) D (P 2 ~ Af?)((p - Kf - Ml) U-=i((P ~ KiY -ml)' 

The coefficients c n ^{e) are rational functions of external momenta and polarization vectors. They are also 
rational functions of e, where D = 4 — 2e by using dimensional regularization. Since these master integrals 
are relatively well understood the problem is reduced to the computation of these coefficients. 

There is another equivalent expression for (|1.1|) that is often used in practical calculations as 



A = c 5ii (e = 0)/ 5 D ' (i) [l] + c 4il (e = 0)/f + c^e = 0)lf 

+ c 2 ,i(e = 0)/ 2 D ' W [l] + c M (e = 0)/f' (i) [l] + (rational part) + 0(e), (1.3) 



1 Here we consider the general massive one loop amplitude, i.e. the propagator is massive. The massless case can be 
obtained by taking Mi = Mi = m,i — 0. Notice that in the massless case the set of master integrals is simpler, since there are 
no one-point functions (tadpoles) or two-point functions with massless external momenta (massless bubbles). 
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where coefficients of master integrals no longer depend on e, but extra contribution - rational part - 
appears. The part expanded by master integrals is usually called cut-constructible part, because coefficients 
Cn,«(e = 0) can be calculated by the pure four-dimensional unitarity cut method 2 |5|, [|. In last few years, 
especially motivated by the twistor string theory 0, techniques for one-loop amplitude calculations are 
much developed and various amplitudes up to the full six-gluon cut-constructible part have been calculated 
|, §, 0, 0, EH, U, © H 0> IH> IH H- However, the rational part is totally lost in the pure four- 
dimensional unitarity method. It has to appeal to some other methods to be computed, either by using the 



improved traditional tensor reduction methods [21, pffi| , or the unitarity bootstrap approach [23, 24, 25], 
or other methods. Recently, an automated package BlackHat has been developed to computer the full 



one-loop amplitude, by combining the four-dimensional unitarity method and the bootstrap approach [26]. 
See |2j], [28| for a full review. 

Another very noticeable development is the method of OPP-reduction |2^] (based on the technique 
in |30]]). Inspired by this work, several very efficient numerical methods have been developed |3l], |32|, 
33, 34, 35, |3q| . In [p7| , a fully automated one-loop N-gluon generator — Rocket has been developed by 
implementing the OPP-reduction and D-dimensional unitarity method. Some analytic techniques were 
also developed in |58L |39[] . The rational part can also be determined in the OPP approach, by keeping the 

9- 



full D-dimensional dependence in all terms [33, 34, 37, 

The .D-dimensional unitarity method is an extension of the four-dimensional unitarity method. This 
idea origins from the fact that a null momentum in (4 — 2e)-dimension can be equivalent to a massive 
momentum in four-dimension [^l], 42, 43]. The power of D-dimensional unitarity method is that it can 
calculate the full amplitude. In another word, it makes the rational parts (also O(e)-terms) cut-constructible 
too. This can be understood by the following expansion of the amplitude 



A = c 5 > 2 ) ® I 5 AW [1] + C4iJ (/i 2 ) ® lf W [l] + c 3 , 4 Gu 2 ) ® 

+ C 2 ,i{f) ® I? + C M GU 2 ) ® /f' (i) [l], 



(1-4) 



where the coefficients c n) j(/i ) can be obtained by using D-dimensional unitarity method [44, 45, 46, 47 



(see also [33, j3f|). The operation "(g)" is defined as 



(1.5) 



The coefficients c ni j(// 2 ) are polynomials of fj, 2 p0|, 49, H3] 3 , where fl is the — 2e component of the (4 — 2e 



2 If using four-dimensional unitarity method, we would not obtain pentagon terms. However, the pentagon and box master 
integrals are not independent: the scalar pentagon can be expressed as a sum of five scalar boxes up to 0(e) corrections Q. 
Therefore, the box terms obtained by four-dimensional unitarity cut include also the pentagon contributions up to O(e). 

3 A complete proof is given in [^9| for massless case, while using this proof, in Appendix of has briefly discussed the 
generalization to massive case. Although appeared in arXiv earlier than [^9), as mentioned in citation of pq|, p^ ] should 
be taken first logically. We emphasize that it's necessary to extract the pentagon terms so that the box coefficients can be 
polynomial of ^ . 
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dimensional internal loop momentum p. By using following relation [50, 42] 



K/* 2 )*] = ^z^^ +2k [l] = ~* m + 0(e), (1.6) 

all /^-dependent terms can be converted to terms of dimensionally shifted master integrals, i.e. I^ +2k [l]. 
So the expansion can be reexpressed as 

A = c M (/x 2 = 0)/ 5 D ' W [l] + c 4ji (M 2 = 0)/ 4 D ' W [l] + c 3 > 2 = 0)/ 3 D ' (i) [l] 

+ C2t i(n 2 = 0)/ 2 D ' (i) [l] + c^ip 2 = 0)/f' (i) [l] + (dimensionally shifted integrals). (1.7) 

The coefficients c Uj i(fj, 2 = 0) are the same as the coefficients c n ^(e = 0) which can be obtained by pure 
four-dimensional unitarity method. The terms of dimensionally shifted integrals will produce the rational 
part and 0(e) terms. The reason is that the coefficients of the dimensionally shifted master integrals are 
always of C(e)-order (which is easy to see from (|Q|)), so only the divergent parts of the integrals could 
give finite contributions, which are the rational parts. 



In the series of work [E6L 47, [49|| , formulas have been given to compute the coefficients c n) j(//. 2 ) 
(except the tadpoles and massless bubbles) 4 . By using these formulas, the pentagon, box, triangle and 
bubble coefficients can be read directly from the tree level input without evaluating any integrals. However, 
a potential weakness for these formulas is that the tree level inputs are required to have no spurious 
pole. This is in contrary with a general observation: while the tree level inputs obtained by BCFW on- 



shell recursion relation 51] are in a very compact form, they usually contain spurious poles. To take full 
advantage of the power of on-shell recursion relation, it is necessary to generalize the formulas so that they 
are directly applicable for the tree level input in the presence of spurious poles. This is the main goal of 
the present paper. While solving this problem, we are able to give very general algebraic expressions for 
most general tree level inputs. 

The paper is organized as follows. In the next subsection, we list our main results. In Section 2, we 
give a brief review of the D-dimensional unitarity method and define various forms of tree level input. 
In Section 3, we study some properties of the formulas from which we can gain some perspective on 
how to make generalization. Then in Section 4, we present the generalized compact formulas, with a 
rigorous proof. Two further problems are solved in Section 5. First we give a general pentagon formula 
by using the quintuple-cut method. Then we simplify the it-dependence of the box formula based on the 
previous result. In Section 6, we implement the new formulas to compute the full five gluon amplitude 
A(l~ , 2 + , 3 + , 4 + , 5 + ). A summary is given in Section 7. In Appendix A, we summarize the formulas in 
previous works. In Appendix B, we show the new pentagon formula is consistent with the previous one. 
Then in Appendix C, an equivalent expression is given for £ij which is used in box formula. Finally, we 
give explicitly the /it 2 terms of the box coefficients in Section 6. 



4 The coefficients of cut-free functions like tadpoles and massless bubbles can't be obtained via double-cut unitarity. In the 
massless case, these coefficients identically vanish. But in massive case, other methods are needed to evaluate them, for example, 
by considering the known divergent behavior of the amplitude ^| , or with alternative techniques [^, EiL |3s], ^] . 
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1.1 Summary of our results 

In this part, we list our main results in this paper. The definitions of various variables and functions can 
be found in later sections. 

The starting point is following double cut integral 



-i(47T 



.D/2 



d D p 



S (+)( p 2 _ Ml) <5«((p - K) 2 - Ml) T( N \7), 



where 



T^ N \i) = A^il) x A% ec (£). 



(1.8) 



11-9) 



The T^ N \£) can be calculated by any method, for example Feynman diagram, off-shell recursion relation 
p2| or BCFW method J51] ], and allow the presence of spurious poles. N is the degree for p and will be 
defined in Section 2. In our notation, p = £ + fl and 

-0y/T=Z ( P x ~ x 



(e\K\£] 



AA 



K 2 



K ) — a- 



K ■ Rr 



AA 



K 2 



K 



(1.10) 



where P xx = \£) \£] and u,a,(3 are given by ( 2.14 ) and ( 2.16| ). We can also rewrite = —2K ■ P xx . 

We should treat \£) and \£] as independent variables when we make replacements given below. 
Now we list our results (Figure 1 shows the graphs that correspond to various formulas): 



• (a) Pentagon: It is given by 

V^K^K^K^K] = rW(% r) ) • Di (I^j^Dj (% j, r ) ) D r (% )i)r ) ) . 
where £(ij t r) is given by (5A) and Di(£) is given by ( [2.22| ). 

• (b) Box: It is given by 

BoAK h K h K] = \ frW(4) • A(4)^(4) " E PeQ[Ki > K i> K r> K] ) 

1 \ r#j Dr{£ij) J 

+{Pji,l Pji,2} 

where £ij are given by |U1|) or @, while \P ji<a ) and \Pj ] are given by (|A.13| ) 

• (c) Triangle: The triangle coefficient is 



(1.11) 



\t] - \Pji,t\ 
W - \Pn,i) 



(1.12) 



Tvi[K s ,K] = - 



(K 



2\N+1 



2 ^f 3 ^T^)N+l^_ Aq 2 K 2 ) N+l {N + !), { p sA Pg2 

^( m\i) N+1 TiN){I) . D(I) 

dr N+1 \ (K 2 ) N+1 K ' K ' 



N+l 



\£\ - \Qs(u)\£) +{P S ,1^P S ,2} 
\£) -> |P s ,i - rP 3 , 2 > , 



(1.13) 
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p — — 



P-K, 



p - 



/p - K r 



Pen[Ki, Kj, K r , K] 



BokIK^K^K] 




p- K 



TA[K it K\ 




Figure 1: The figures corresponding to various formulas are shown. Notice the order of the propagators can be 
changed. Our convention is that all external momenta are taken outgoing. 



where £ is given by ( |1.10|) , Q s (u) is given by ( A. 2 ), and P S)Q is given by ( A. 14 ). 

It is worth to emphasize the ordering of replacement. We need to replace \£] — > \Q s (u) \£) first where 
\£) shows up. Then we replace \£) — > \P s ,i — tP s $)- To denote this ordering, we have put the first 
replacement above the second one. Similar understanding should apply for later formulas. 

(d) Bubble:The bubble coefficient is 5 

r=l a=q 



N 

Bub[K] = {K 2 ) N+l 



k N 



d 



N 



dr N 



q=0 

(2ry • K) t+1 



q\ dsi 



(t + 1)(K*)W N\[ V \rjK\r,}" (£ V ) N+1 { (# 2 )" { >) 



(r,a-q,2), 
N,N—a 1 



(1.14) 



\£] - \K + s V \£) 



1.15) 



d b 



b\ {(3^T—u)b+l ^- Aq 2 K 2 )b +l (p fA p r jb dT b 



P 



it+1 



t + l^|K|P ril ]' +i (%-ftT^) 



M -» |X + sr,K> 

K> -> \Pr,l - rP r , 2 ) 



(1.16) 



1 



6! (/3VT^) b+1 (V"4g^ 2 ) fe+1 (P r ,l Pr,2) 6 ^ 6 
iJV+1 N 



1 (%IP, 



r,2 



(€|Q r («)r7K) f 



t + 1 (£\K\P r , 2 ] t+i (£\r)K\£) 



,JV+1 



JV+l 



-TW(£) ■ D r (£) 



K) -> |Pr,a - TP r ,l> 



(1.17) 



where £ is given by (p-lOp , Q r (u) is given by ( |A.2j ), P rja is given by ( |A.15| ), and 77,77 are arbitrary, 
generically chosen null vectors. The summation over r includes all the D r {£) that appear in the 
denominator of T{£). 



J Here we give the formula with general choice of auxiliary null momentum r\. For some special choice, the formula may be 
modified, which can be found at the end of Section 4. 
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If N < —2, there are contributions only from boxes and pentagons. If N > —1, contributions from 
triangles will kick in, and finally if N > 0, bubble contributions show up as well. Similar observation has 
been made in [|[. 

Our formula should be suitable for both analytic and numerical calculations. For numerical calculation, 
since our result should be polynomial of u, so we should have a way to get coefficients of u a numerically. 
This can be done by using discrete Fourier transformation method ||35|| . 



2. D-dimensional unitarity method 

In this section, we briefly review the .D-dimensional unitarity method [^] and define various forms of tree 
level input. In the process, we also establish our conventions and give definitions of some variables and 
functions. Some preliminary knowledge about spinor formalism [53] and color decompositions [54] can be 



found in two classic reviews [55, |56[. Here we only emphasize that, we use the "twistor" convention for the 



square bracket [i j], so that 2ki ■ kj = (i j) [i j}. 

Let's start with the following expression for a double-cut (iT 2 -channel) integral: 

^] ^-channel = ~ W ^ / " M?) 5((p - K? - M|) T(p). (2.1) 



Compared with the full amplitude expressed in (1.4), i.e. 



A = J2 c n ^ 2 )®I°^[ll (2.2) 

n,i 

the cut integral Q2.1| ) can be understood as 

n,i 

where the cut master integral is 

r r^rill I f d D p 5(p 2 -Mf) 5((p-Kf-Mj) 

[n[JJI ^ 2 - c "" 1 j 7(2^ HEfUp-Kty-m*) ' [ ] 

Those master integral, which don't contain two propagators that correspond to the two 5-functions, will 
not appear in the cut integral fl2,3| ). By evaluating different cut integrals, we can get all coefficients c ni «(^ 2 ). 

In next subsection, we will briefly review how to simplify double cut phase-space integration. This 
technique is the foundation of our method. After these transformations, we can describe the structure of 
cut integrands, from which it is possible to read off coefficients directly. 
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2.1 Simplify phase-space integral 

The internal loop momentum p is (D = 4 — 2e)-dimensional and can be decomposed as 

p = £ + ft J d 4 ~ 2e p = J <T 2e iJL j d 4 £ , (2.5) 

where £ is 4-dimensional and ft is (— 2e)-dimensional. We can further decompose £ as 

£ = £ + zK, £ 2 = 0; J d 4 £ = J dz d 4 £ 5( +) (£ 2 ) (2£- K), (2.6) 

where K is the momentum across the cut. I is a massless 4- momentum, and can be expressed with spinor 
variables as J57|] 

I = tP xJ , P xl = XX = \£) \£] ; J d 4 £ 5 {+) (£ 2 ) = J (£ d£) [£ d£] J t dt. (2.7) 

The corresponding measure has also been given. 

The Lorentz-invariant phase-space (LIPS) of a double cut is defined by inserting two (^-functions 
representing the cut conditions: 

d 4-2 e$ = J d 4-2e p § ^ p 2 _ M 2 )5((p - K) 2 - Ml) (2.8) 
(47r)£ f dfi 2 (fi 2 )~ l ~ e f d 4 £8(£ 2 -fi 2 - M 2 ) 8((£- K) 2 - fi 2 - M 2 ) 



T(-e) 

= ^1 J dfl * (V 2 )- 1 - 6 J dz d 4 £5^(£ 2 ) (21 -K) 

x5(z{l - z)K 2 + z(M\ - Ml) - Ml - fx 2 ) 5(-2£ ■ K + (1 - 2z)K 2 + M x 2 - M 2 2 ). 
We can firstly perform the integral over z with 5-function S(z(l — z)K 2 + z(Mf — Ml) — Mf — fi 2 ) to reach 

J d 4 ~ 2e <f> = ^1 J dfi 2 (p 2 )- 1 -' J d 4 £ (£ 2 ) 5(-2£ -K + (l- 2z)K 2 + M 2 - Ml) (2.9) 
J dfi 2 (fi 2 )- 1 ' 6 J (£ d£) [£ d£] j t dt 5(-2tK ■ P xx + (1 - 2z)K 2 + Mf - Ml), 



(47T) 



r(-c) 

where z is fixed by solving the <5-function as 



_ (K 2 + Mf-M*)- y /A[K,M 1 ,M 2 ]-<lK*v? 

Z ~ 2K* ' [ ' 

with 

A[K, Mi, M 2 ] = (K 2 ) 2 + (Mf) 2 + (Ml) 2 - 2K 2 M\ - 2K 2 M 2 2 - 2M 1 2 M 2 2 . (2.11) 
Then we can perform the integral over t with the remaining (5-function to reach 

, (2 , 2) 
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and t is 



(1 - 2z)K 2 + Ml - Ml 



2K ■ p xx 



(2.13) 



For convenience, one can redefine the // 2 -integral measure as 



2/ 2\-l-e 



d/J, (n 



f A[K,M 1 ,M 2 ] 



-e /•! 



(in it" 



-l-e 



where the relation between u and /i 2 is given by 



4*V 



A[K,M u M 2 y 



(2.14) 



Then we can rewrite z, t as 



K 2 



2K ' p x~x 



(2.15) 



where 



a 



i^ 2 + M x 2 - Mf 
K 2 







VA[JT,Mi,M 2 ] 
if 2 



(2.16) 



Notice that when Mi = M 2 = we have a = /3 = 1, thus reproducing the massless case. A useful relation 
between z and u is the following: 



. M 2 - M? „ , 

( X " 2z ) + = ^v 7 ! 3 ^. 



(2.17) 



After these evaluations, the cut integral (2.1) is transformed to 
(4vr) e 



C[A] 



m D l 2 T{-e) 
(Any 



?7T 



D/2Y(- 



d^ 2 (u 2 )- 1 "^ 

A[K,M 1 ,M 2 } ^-' 
4K 2 



(^)[^] (1 " 2 lf t M r M2 V(p) 



(2* • r,) 2 



K 2 



(2K ■ P x ~ x ) 2 



Here T(p) should be interpreted as 



T{p) = T{1 , u 2 ) = T(ti> r + zK , u 2 ) = T(|*> , \t] , u 2 



T(p). (2.18) 



(2.19) 



with z and t defined as in d^JEj) . In the following, we may write Tip) as T(£) where the dependence on 
u 2 is implicitly included, with i being interpreted as 



K 2 

tp ~ + Z K = 



-pVT^u ( P xJ tt^K ] - a „ n xx K 



K 2 



K 2 



(2.20) 
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2.2 Standard form of the input 

For a double cut of physical amplitudes, it is always possible (for example from Feynman rule) to write 
T(p) as a sum of terms of the following form 

T (n) ,r* _ ll7=i fc (- 2 l- Pj) , . 

where Di{l) is the propagator 

Di(I) = (I- Kif - /i 2 - m? = -2? ■ Ki + Kf + Ml - mf. (2.22) 

We emphasize that the relation £ 2 = M 2 + fj? from the cut constraint is always assumed in the input. 

We will call the form of T(p) like ( |2.21| ) as the standard form, where there are only Di(£) in the 
denominator. We will frequently use an important quantity - the degree for p (or equivalently for €). The 
degree of T{p) is defined as the degree of numerator minus the degree of denominator. For the standard 
form fl4.9|) , the degree is n, which was shown explicitly in the superscript. 

After simplifying the phase-space integration, the cut integral for the standard form can be written in 
the form of 



iir 1 - 



V j j Q 7 v u\m + nti a\Qi\i] 



r D/2 r(-e) V 4^ J J J (£\K\£] n+2 Ui=iWQM 

where 



Qi(u) = -/3(V1 - u) ( -^i ^2"^ J ~ ( a ~ ^2 ^2 ) ^ 



and a, (5 are given by ( [2.16 ). We should notice the equivalent relations at the integrand level as 



- 2? ' p ' = W • A(?) = WW) mM • (2 ' 24) 

and therefore 

For the cut integral with standard form input fl2.23| ), we can read off coefficients directly by using 
several formulas, which are reviewed in the Appendix [A|. 
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2.3 General form of tree level input and spurious pole 

We now consider the general form of inputs. In unitarity method, the input needed is a collection of 
(on-shell) tree level amplitudes. By using the on-shell recursion relations, very compact expressions for the 
tree level input can be obtained. However, such compact expressions usually contain some spurious poles. 

A general form of T(p) that with a spurious pole can be written as (for simplicity, we only consider 
the case of single spurious pole, multi-spurious pole is a trivial generalization) 



(2.26) 



where q are coefficients that do not depend on £, but may depend on fi 2 . Sd is the spurious pole with 
degree d and can be generally expressed as 



S d = s + s i(~ 2i ■ Vi,i) + Yl »<i.*a(-2^ • V 2M )(-2l ■ V 2)i2 ) + ... 

i ii,i 2 

+ £ s ll ,..., ld {-2l-Vd, ll ){-2l-Vd,i 2 )...{-2l-V d , ld ). 



The degree of Tgp (£) is N, which is defined as the maximum of nt 

N = Max{nJ. 

After putting into ( p. 23 ), Tsp \l) can be rewritten by using the relation fl2.24|) as: 

7$°® = E^- {K2)nt UT = lh+d m]i] 



where 



Sd = so 



(£\K\e] d ^ {e\K\£] d - s 



(K 2 ) d 

+ Y s h. 



\K\£] 



d-2 



d-2 ° l ^2 



i\V dth \£ (£\V d>i2 \£ ...(£\V d ^ d \£ ] , 



(2.27) 



(2.28) 



(2.29) 



(2.30) 



and Vt.i z is obtained from Vt,i z just as that Rj obtained from Pj. 

For a physical amplitude, poles other than Di{£) are all spurious. It is always possible to regroup 
terms into an expression free of spurious poles, i.e. a sum of standard forms: 



T-rn r +fc 



^■Pj) _ v {K 2 r ngi^wj 



(2.31) 



b r are coefficients that do not depend on £, but may depend on // 2 . The degree, now the maximum of n r , 
should also be N. 
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An important observation is that: to go from Tgp\l) to Tn^p(t), we only need to apply some spinor 
identities to remove the spurious pole, while I must satisfy the two on-shell conditions: 

I 2 -l? - M\ = 0, (7- K) 2 -ii 2 - Ml = 0. (2.32) 

This observation will be important in the proof of the generalization in Section ||. 

3. Looking at the formulas 

The formulas we have known (collected in Appendix |A|) are applicable for the standard form T s f(p), in 
which the key requirement is that there is no spurious pole, i.e. only propagators D (I) appear in the 
denominator. One may expect to remove the spurious pole in a general input to reach the standard form, 
since it's possible in principle. However, in general practice, this procedure is very complicated and we 
do not have well-defined algorithm to do that. Besides, the most compact expression of the tree level 
input usually has spurious poles. After removing the spurious poles, the expression may expand quite a 
lot, which leads the computation to be less efficient. Therefore, it is important to be able to generalize 
formulas to the case with spurious pole, while keeping the compactness of the tree input. 

In the following, we will study the previous known formulas for the standard forms input from a new 
point of view, by which we can gain some insights on how to generalize the formulas to the case with 
spurious poles. 

3.1 Splitting of the integrand 

We start from the integrand in (|2,23j ) 



7 _ (K 2 r+ i n^mm 
~ {i\K\i\ n+2 nf =1 ' 



(3.1) 



Then we do the splitting to reach 

k 



I ~ 1 

1 = I> (A Wi <wi + S Gi(A ' A) (wP 7 ' (3 - 2) 



=1 * ' 3=0 



where 



(£\KQ l \£) n+1 nt Wi (£\QtQi 



t^j^j - (^f^^^L , (3-4) 



j=0 \-|"P-J q =o 
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with 



B n ,t( s ) 



(£\K\£] 2+t (£\rjK\£) n ft* (£\Q P (K + srj)\t) 



(3.5) 



while the Gj term will 



The box and triangle contributions will come from the F% term in ( p.2| 
contribute to bubble. Taking the integrand into the cut integral, and extracting the residues of various 
poles, we can get the final coefficients. With the above splitting, we will know all coefficients by putting 
them into following expressions: 

2 W x 



C[Q u Qj,K\ 



+ {p jifl «-► p^} 



(3.6) 



for boxes defined by Qi, Qj, K, 
{K 2 ) n+l 



C[Q S ,K] 



2(VAI)™ +1 (n + l)\(P s>1 P S!2 



n+l 



d 



'n+l 



n?=r ( p s,i - TP s ,2\RjQs\p s ,i - tp ; 



8,2} 



dT n+l 



,nti, ¥s (Ps,i-rP s , 2 \Q t Q s \P s , 1 
for triangles defined by Q S ,K and finally for bubble 

(-!)" ,1'' 



s,27 



(3.7) 



C[K] = (K 



2\n+l 



E 

q=0 



dsi 



B 



(0) 



n,n—q 



w+EE 



(r;o-g;l) 



g (r;a- 9 ;2) 




(3.8) 



r=l a=ij 

where definition of various variables and functions can be found in Appendix |A]. The key observation is 
that all coefficients are just the sum or difference of relative residues of pole. In Fi(\), the poles 
from (£\QtQi\£) contribute to box, and those from (£\KQi\£) contribute to triangle. Bubbles are from the 



poles in Gj, or more precisely, B^\(s) is from the poles (£\r]K\£}, and B^f' x> (s), B^'"'^'(s) are from the 
poles (£\Q P (K + srj)\l). Details of the derivations can be found in Jiq , 47 1 

Let us recall how do we get above formulas. Starting from standard input, we algebraically simplify it 
by using splitting technique. Then we perform phase space integrations carefully and add all contributions 
from physical poles together. If now the input has spurious poles, the same method tells us that we need 
to perform phase space integration with spurious poles. However, we can not do that in general. The 
reason is following: (1) The form of spurious poles can be arbitrary and we do not know how to write them 
down; (2) To perform the phase space integration, we need to write it into total derivative and then take 
residues. For a given spurious pole, we do not know how to perform both manipulations. 

Because these reasons, it seems impossible to find results by standard method. However, considering 
these poles are all spurious, this means that even we are able to extract the residues of spurious poles 
term by term, their sum should vanish from physical point of view. In this sense, we can expect that the 
formulas can be generalized to the general form without many modifications. Then the question is how to 
do it? Which new method will enable us to do it? 



(r,b,l). 
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To be able to proceed, we find that we need to take a new perspective in understanding above results. 
Remembering that when the spurious poles show up, there are at least two terms with spurious pole, so 
the added result will cancel out. In another word, we will have following equations (we assume there is 
only one term of standard form for simplicity): 

T{p) = -^r (3-9) 

nti aw 



where is given in (p7|)and {ij means the set {i±,...,it}. Among these three forms, the first form is 



the one without the spurious pole, while the third form is the one obtained by using on-shell tree level 
recursion relation. The second one is the bridge bring the first one to the third one. 
After putting into ( |2.23| ), we have following three forms for the integrand 

r K 2)n+i Y[^(t\RM 

(i\m n+2 nti mm 

{K 2 T+ i+t- d n f z= i (#m,ii] u]=i mm 

^ Si i-*2,-,ii , ]n+2+ t-d c n e \ rrfc io\r> \ff\ (3.13) 



{K 2 )nt+ i- d uT=i kt mm (3U) 

" V Ct (*\ K\i] nt+2 ~ d s d [\t) , \£]] nfei (MM 

which should be considered as function of \£) and \£]. It is important that we have treated \£) and 
\£] as independent variables. 

3.2 The splitting result for box and triangle 

The next step is to find the expression contributing to various coefficients. 

Let us start from the box and triangle. From the form ( |3.1| ) we can see that the relative expression 
(i.e., for box and triangle having Qi) is 



(I) = ^ ) ll J= i v^^i 

{£\KQi\i) n+x I[t=i,t^\QtQi\£) 

Comparing ( |3.15 ) and (|3.1| ) we found that it is obtained from (3T) by following manipulations: 



(a) Multiplying (£\K\£j and {£\Qi\£j at Q; 
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• (b) Then replacing \£] — * Qi \£). 

Using this observation, we can get two results from ( |3.13 ) and ( |3.14 ) as following 

{K 2 )n+ i + t^ Ul=x (t\v t> uQi\l) U]=i (i\RjQi\t) 



F; 



(II) 



s h,i2,-- 
t&t} 



(t\KQi\l) 



n+l+t-d 



(3.16) 



and 



F, 



(III) 



(K 



2\n t +l-d 



(e\KQi\Q 



nt+l—d 



s d [\t),Qi\£)]Uti (t\QiQi\t 



(3.17) 



It is worth to give a remark regarding the expression (3.17). In this formula we have factor (l\QiQi\i) = 
in numerator. If there is same factor {£\QiQi\£) in denominator, then both factors will cancel each other, 
but if there is no such factor in denominator, we know immediately that this term contributes zero. A 



subtle point for box formula will be discussed in Section 3.4 



Since three expressions ( |3.12| ), ( |3.13| ) and ( 3.14| ) are equal to each other at the algebraic level. After 
the same algebraic manipulation, we should have that three expressions ( PHI ), fl3~l6D and (gl7|) are equal 
to each other too. 

Now we can read out the expression for the box and triangle coefficients. 

3.3 Read the coefficients 

The coefficients of box from ( |3.15| ) is 



1 (K 



2\n+2 



2 (£\KQ t \e) n+2 U k t=1 . 



+ 



1 (K 



2\n+2 



\i)^\pn,i) 



2 {£\KQ i \£) n+2 Y\ k t ^ i .{£\Q t Q i 

which is equivalent to (|3.6| ) by noticing that 

(P ]hl \XQi(u)\P jiA ) = (P jiA \X\P jit2 } 
(Pji^X'Qiiu)]^) {Pji,i\X'\Pjid ' 

Compared to (3.15), the coefficient can be obtained by following algebraic actions: 

• 00 Multiplying tffgf ; 

• (b) Sum up two terms with \£) — > and \£) — > \Pji,2(u)). 



(3.18) 



(3.19) 
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Now it is obviously we should apply same actions to (gig) and fljUTD to get 



C[Qi,Qj,K]W 



E g n,»2,- 
9 



2\n+2+t-c! 



n+2+t-d 



+ 



KHI^i,l> 



5 d [i€),QiW]nti,7^^w^ 



(3.20) 



and 



C[Q h Qj,K] 



(in) 



2\n t +2-d 



, 2 ^|iTQ i |£) nt+2 - d 



U n s L + i kt (t\RsQi\i) (i\QiQi\i) (t\QjQi\e) 



U n s L\ kt (£\R s Qi\i) (£\QiQi\i) {i\QjQi 



+ 



(3.21) 



Formula ( 3.21 ) is the algebraic expression for box in the presence of spurious poles. We got it by various 
algebraic replacements starting from input tree-level amplitudes. We can do similar algebraic replacements 
to get algebraic expressions for triangles and bubbles too. For later two cases, there are two complicities 
compared to box. First there is operation of taking derivatives. Second, formulas depend explicitly on 
the degree of the tree level inputs. These two issues can be solved with some considerations. In the next 
section, we will give a rigorous proof for the generalization of these formulas. The final formulas will be 
equivalent to those obtained by naive substitutions. Before going to that, let's have a look at an example 
demonstrating some subtle point regarding the box coefficients. 



3.4 Subtle point regarding to box coefficients 

In the standard form we know that if one term contributes to triangle coefficient C[Qi, K], it must have Di 
in its denominator. Similarly if one term contributes to box coefficients C[Qi, Qj,K], it must have Di, Dj 
in its denominator. 

If the form is not standard, i.e., having spurious pole, above observation will be different. For triangle, 
above observation is not modified, i.e., if one term contributes to triangle coefficient C[Qi, K], it must have 
D{(€) in its denominator. However, for the box, the above observation is not true anymore: one term can 
contribute to C[Qi, Qj, K] even it has only Di or Dj in its denominator (but at least one of two). 

Let us consider one simple example. It is given by 
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for cut K\2- Here we do everything in pure 4D so we do not have u floating around. It is obviously the 
only contribution is box and no others. However we can rewrite it as 

r 



{21 ■ K 41 ){£ - KrY {2£ ■ K 4X ){£ + K 4 ) 2 ' 

71 = " ' {-2£ ■ K 41 ){£ - 72 = {-2£-K 4l ) {£ + K 4 ) 2 ' (3 ' 23) 

where {21 ■ K 4 \) is the spurious pole. Now let us apply our general method to find various coefficients. 
Before doing so, let us list following quantities (we have set z = u = 0) 



Qi = -K u Q 2 = K 4 , V = -K. 



41 • 



First let us start from box. Spurious pole shows up in two terms, one with propag ator {£-K x ) 2 and 
another one with propagator {£ + K 4 ) 2 , thus we need to consider all possible combinations of boxes 
in these two sets. In our simple case, these is only one option: box with {£ — Ki) 2 and {£ + K 4 ) 2 . There 
are two terms, both with kt = 1, d = 1, and i%t = —2. 

Let us apply formula ( 3.21| ) with the choice i = l,j = 2. For the second term the contribution is 





i\VQx\t) J z \{i\VQi, , 

This can be seen immediately since the denominator does not have D\ propagator. For the second term 
71, the contribution is 




Adding two contributions together we see that indeed we have reproduced the right box coefficient. 

However, there is an subtle point regarding the calculation of T\ in this example. In fact, when we put 
in the solution P±,P2 into (£\Q2Qi\£) and (^£\VQi\i^ , both are zero, so we need to take proper limit. The 

point is that even £ be arbitrary we have {£\Q2Qi\£) = (£\VQi\£S and (£\QiQi\£) = 0. 

Now let us move to triangle, since n = —2 we see that triangle coefficients are zero as we familiar with 
(when the power of derivative is negative, we should take it as inserting r' a ' — ► 0). For bubble, we have 
n = —2, again the contribution is zero. 

It is worth to emphasize again that above subtle point is related to box coefficients only. 



4. Formulas with spurious pole 

Now we generalize the formulas with a rigorous proof. As a byproduct, we will obtain more compact 
expressions for the formulas. We will take two steps. First we reformulate the formulas for the standard 
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form input into a more compact form, in which the tree level input is inserted directly. Then with this new 
presentation, it will be easy to generalize the formulas to the tree level input in the presence of spurious 
poles. 



4.1 Reformulate the formulas 

We first reformulate the formulas for the tree level input being standard form. Using the relation ( p. 24 ) 



-21 ■ P, 



K 2 



j 



{i\K\£\ 

we can reformulate the box formula (A. 5) as 



{e\Rj\e\, Di(e) 



K'- 



[£\K\£] 



no ok]- 1 ( ^ 2 ) w+2 rfc w W 

C[Q U Q V K] - - I ——^— k 



2 \(i\K\e] n+2 nUt^(t\Qt\t] 



\t] - \Pji,2] + {Pji,l <-> Pji.2] 

\l) - \PJi,l) 



\i] -> \Pji,2] + i P jiA ^ P ji,l} 

\i) - \PJi,l) 



Similarly, for the triangle formula (A. 7), we can get 
(K 2 ) n+1 1 



C[Q a ,K] 



2(VA^ n+1 (n + l)!(P,,i P s , 2 ) n+1 

■DM 



( {£\m i+1 T (n) , 



dT n+l I (^2) 



n+1 sf 



W - 10. K> + {P s ,i *-> P s ,2} 

\£) -> \Ps,l - rP s , 2 > , 



and for bubble (P)-(A32), we have 



«£](-) = ^ 



(2rj • K) t+l 



\t] -» \K + s V \i) 



(4.1) 



(4.2) 



(4.3) 



(-i) 



6+1 



WCV^)^ 1 (Pr,l Pr,2) b dT b 



1 ( (i\K\e\ n+1 n-W 



t + l(£\K\P rA ] t+1 (£\ V K\£) 



\n+l 



|/| _ |X + s^|£> 
K> - l-Pr.l ~ rP r , 2 ) 



(4.4) 



(-1) 



6+1 



6!(VA;)^+i (P rA P r o) b dr b 



1 <%l^r,: 



ii+1 



(£\Q rV \l) b {£\K\£T +l {n) ~ ~ 

t+l IB\„v\e\n+\ \ (K2\n+1 2 sf V-l'^rV^) 



t + l(£\K\P r , 2 ] t+1 (£\ V K\£) n+1 \ (K 2 ) n 



\i] |K + sr^> 

\£) ~* \Pr,2 - rP,,l> 



(4.5) 



- 18 - 



We can find that the formulas actually preserve the structure of tree level input {€) quite well. 
The main structure of the formulas is that: first make a substitution for \£], then for \£) , and finally do 
some operations at algebraic level. 

The substitution for \£) and \£] is equivalent to a substitution for £, by using the relation fl2.20| ) 

K 2 

_ .1 —Hi i 

xx 



[£\K\£] 



K 2 



a 



XX K 



K 2 



(4.6) 



It is not difficult to find that, with these substitutions, £ satisfies the two on-shell conditions ( [2,32 ), because 
P^r under the substitution satisfies the massless condition. For box this is obvious. For triangle and bubble, 
we have 



P 2 X ~ oc (£\£) = (P s>m - rP S)2(1) |P s , 1(2) - rP S)2(1) ) = -r((P S)1 P s , 2 ) + (P S)2 P Sjl )) = 0. 
4.2 Generalize to general tree level input 



(4.7) 



To generalize the formulas, a naive conjecture is that: for a general form of T{p) with degree N, the 
formulas Q4.1|)-(fOD are unchanged, but only with the following substitution 



N, 



where 



(K 2 r Y\7=i t+d mw) 



(4.8) 



or equivalently in a form without spurious pole 



nf=iAW 



The degree N of T(p) is defined as the maximum of rit or n r : 

N = Max{n 4 } = Max{n r }. 



As mentioned at the end of Section 2J3, the important observation is that: to go from Tsp\l) to Tnsp {£), 
we only need to apply some spinor identities to remove the spurious pole, while £ must satisfy the two 
on-shell conditions ( 2.32 ). 

We now want to prove the conjecture. First we can see that if the formulas is true for the substitution 



(4.9) 



(4.10) 



then it will also be true for the substitution 



(4.11) 



(4.12) 
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This is because in the formulas the substitutions for £ do not break the on-shell conditions fl2,32|) as shown 
in the end of last subsection, and also by the following simple observation: if two functions fi(p) and f2(p) 
are equivalent by some algebraic operations, then it would also be true that 



d n 

-^{h(p)\p-*q{ S )] 



-T-[h(jp)\ P ^ q (s)] 



(4.13) 



+ o ds n ' 

Thus to prove the conjecture, we only need to prove for the case ([4.11| ). Since T^p \i) is a sum of the 
standard forms with different degrees, and iV is the highest degree, what we need to prove is actually that 
the formulas with lower degree can be reexpressed with a higher degree. We will prove this case by case. 

4.2.1 Formula for box 

It is trivially true for box formula, because the structure of the formula doesn't depend on the degree of 



input. So we get the formula directly by using (4.1) as 



C[Qi,Q 3 ,K] = - {TW-DiW-DM 



\t) 



\P*,i) 



(4.14) 



for a general form of T(t). 

The formula ( 4.14| ) can be understood from another point of view, as the box coefficient obtained by 
using the generalized unitarity method of quadruple cut, where the substitutions for I just correspond to 
the two solutions that are solved from the constraints of quadruple cut. 

4.2.2 Formula for triangle 

For triangle, it is not so trivial as box, because there are operations of taking derivatives that depend on 



the degree of input. From (4.11), we need to prove that 
(K 2 ) N+1 1 



C[Q S ,K] 



2(VAiyv+i (N + l)\(P Sjl P S)2 ) N+1 



\Qs\t) +{Ps,l^Ps,2\ 
\P.,1 ~ rP s , 2 ) , 



E 6 



2(VAl)'v+i („ r + i)i ( Ps)1 p 9j2 )"<-+i 
<1"-+ 1 ( (l\K\tf r+l 



w A m r (nr) (T) d (I) { 

+1 I (#2)n r +l J sf W J 



\t] - 10* K> + {P a ,i <-» P,?} 

\£) -* [P.,l - rP 3 , 2 > 



(4.15) 



We first do the following calculation: 



\K\£] 



\Qs \i) 

\Ps,l - 



r Ps,2) 



t (P Stl P s , 2 ) VA S 



(4.16) 
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Then by using the relation that 



d(r m f(r)) 



dr n 



n\ d(/(r)) 



r^o (n — m)\ dr r ' 



T->0 



(4.17) 



we can find that 



2\n+l 



d n+1 ( (t\K\iy 



2(VAI)" +1 (n + 1)! (P s ,i P s ,2) n+1 dT»+i \ (K*)' 

(K 2 ) n - m+1 1 
2(V5I) ra - m+1 (n - m + 1)! (P Sjl P s , 2 > 

^rt— m+l 



Ifl - 10. W +{Ps,i 
\e) - \p s ,i-tP s , 2 ) 



Ps,2} 



n—m+1 



W - 10. \t) + {P a ,l Ps,2] 

K> - |P S ,1 - rP s , 2 > , 



where the function f(P x j) should be general but without factor |^] in its denominator. So by substi- 
tuting Tnfp{l) in the formula and changing n to N, we have 



C[Q S ,K] 



(K 



2\N+1 



d N+i 



2(^)^+1 (AT + 1)! (P S)1 P a , 2 >" +1 dr"+i 



(#») 



iV+1 



(if 



2\AT+1 



(£\K\£] 
1 



n r +l 



10, \t) + {Ps,l " Ps,2} 

\Ps,l - tP s , 2 ) , 



2(VAiyv+i (N + l)\{P sA P s , 2 ) N+1 

d N+i / v , (4^]^: n"=i fcr (iiwi 



n-=i mom 



(k 



2^n r +l 



2(VAl)^+i (n r + l)!(P M P s , 2 ) 



firs'* m\M\Q.\4 



d T n r +l 



W - 10, K> + {P S) i P s , 2 } 



Ifl - \Q.\Q +{Ps,l Ps,2} 

\i) -> |P Sjl -tP s , 2 ) , 



(4.18) 



which is just what we want to prove. Notice that if one term in T(£) doesn't has D s (£) in the denominator, 
its contribution would be zero since 



[i\Q s \t] 



\Ps,l - rP s ,2> 



. 



(4.19) 
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Therefore we have 



C[Q S ,K] 



2\N+1 



2(VSI) iV+1 (N+l)\(P s ,i P s , 2 



N+l 



dT N + l I {K 2 )N+ 1 1 W 



\t) -> \Ps,l -rP s , 2 > , 



(4.20) 



for general T^ N \p). 

4.2.3 Formulas for bubble 



For bubble, the proof is similar to the case of triangle. Consider the bubble formulas 

^ q\ dsi , 

<j=0 \ r=l a=g 

?(°) ^ R(' r ' a - < ?. 1 )/'„\ „„J K?(r,a-g,2), 



C[tf] = (if 2 ) n+1 



O) + E E K-a'^ 00 - <n a -"a 9 ' 2 ^) 



(4.21) 



where B„„_ g (s), B n ' n _^' (s), and B n ' n _^' (s) are given by ([4.3|) , (4.4) and ([4.51) . We want to generalize 
the formulas to the input (£) . 

(0), 



We first consider B n ' t (s): 



d" 



dr r - 



(2rj-K) 



t+i 



(t + l)(^ 2 )*+ 1 n![7 ? |^|7 ? ]™^7 / ) n+1 V {K 2 ) 



(WW-®® 



As in the case of triangle, we need the calculation 



(£\K\£) 



\i[ - \k + «i\£) = -ts [v\rjK v ] (£ rj) L . , 

\i) -» \K-TTj\v] 



\t) -* \K + s V \i) 
\f) -> \K- rfj\v] 



. (4.22) 



Then by using twice the relation (4.17), we have 



(-1)1 d q f d n 



q\ dsi [dr n 



(2 V -K) 



n — q+l 



(n-q+l){K 2 ) n -i+ l 
f(£\K\£] 



n\[r)\r}K\ri} n (£ rj) n+1 V (K 2 ) 



\£] -> \K + s V \e) 
\ t ) _ |jf - r^|r,] 



2\n- 



E 

9=0 



d T n-m 



(2rj ■ K) 



n— 9+1 



(n-g+ l)(if 2 ) n -« +1 



1 



(K 



2\n— m+1 



E 

<?=0 



(n - m)![r/|^|r/] ?1 - m (£ rj) 



(-1)1 di f 



n— m+1 



10 |Jf + «7|/> 
|€> -+ |X-r5y|7?] 



(2r? • if)' 



ri— m— q+l 



(n — m^r^-fT I^]™ m (£ rj) 



(n-m-q+ \)(K 2 ) n - m -i +1 

n— m+1 



\l] -» |X + sr;|£> 
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So if we make the substitute ( 4.11 ) for the formula of B^\{s) and put into ( (4.21|) , we have 



ff/¥+1 f(-l) j g f <*" 
1 J ^ g! ds" 1 dr^ 

1 



(27? • i^)^"^ 1 



(N-q + l)(K 2 ) N ~i+ 1 



( (i\k\£] n v (^ 2 r- ngtNwg ' 



K> -* |JT - r^r?] 



A? 



r ^2WV+l V l-l) g g f ^ 

1 ^ ^ g! ds« 1 dr^ 

q=0 * v 



r E 



(27? -K) 



(AT-g+l)^ 2 )^^ 1 



\t] -> \K + s V \e) 



E ( K2 ) nr+1 E 



-!)« d 9 f d n 



q=0 



1 



/! ds« dr 



(2t? • K) 



Ur—q+l 



(n r -q+l)(K 2 ) n --i +1 



n r ![7?|r?^|r?]^^7?)"" +1 I nfei 



K) -> 1^ - 



which is just what we need to prove for the generalization. 

With the same procedure, we can also prove it for B^f' 1 \s) and B^f' 2 \s), by using the calculation 



(£\K\£] 



\f] -> \K + sr]\t} 



-s^r?^) 



K>^|Pr,2-rP r ,l) 



Notice that there is a summation over r in ( 4.21] ), i.e. sum over all the D r that appear in the denominator 
of tree level input. We also need to show that if the input has no D r in denominator, it vanish, i.e. 



" ' -1)" d q " 



q\ dsl^ n ' n 

q=0 a=q 



(r,a- 9 ,l(2)). 



, 



(4.23) 



if {£) has no D r in denominator. By substituting ( ]4.4|) or ( ]4.5| ) in it, and using the relation that 



(£\Q r {K + s V )\£) 



w ^ \K + sn\e) = -r(P r ,i P r ,2) y/Ar + s {£\Q r q\t) (\£\ - \K + sr,\e) , (4.24) 

\t) -» |P-,2 - rP,,i> 1 \£) -> |P r>2 - rP rjl > 



we can find this is indeed true, due to the cancellation between the two terms in the right-hand side of 
above relation. 

Therefore, we finally have the bubble formulas 



v 



C[K] = (K 2 ) N+1 £ 



+1 ^(-1)« / (o) 



fc N 



q=0 



q\ dsi 



r=l a=q 



(4.25) 
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-lN 



(2r? • K) 



t+i 



B 



(r,6,l) 



00 



K> — > |-ff— T77|T?] 



(4.26) 
(4.27) 



1 <%l^,i] m 



(r,ft,2) 
7V,t 



i + l(£|^|P r .,i] m (£|r ? ^|^) 7V+1 ^ (i^ 2 )^ 1 
(-l) fe+1 



K) -> |P r ,l - rP r , 2 > 



6!(VA^ 6+1 (n,i^,2) 6 ^ b 



(4.28) 



1 (^P, 



r,2 



lt+1 



t + l(£|^|P r . i2 ] m (£|^|^) 7V+1 I (^ 2 ) 



\£l -» [Jf + «?W 

K> -> |Pr,2 - rP r ,i) 



for general 7"^' (p) . The summation of r is for all the D r (£) that appear in the denominator of T(l) . 



A special choice of r\. As discussed in the Appendix B.3.1 of [47], we can use a special choice of 77: 
choosing rj = K\ in the case where K 2 = 0. By this choice, we have 

1 1 1 



(4.29) 



where a\ is given by (^ 

The general bubble formula for this special choice, i.e. 77 = K\, becomes 6 



N 



C[K] = (K 2 ) N+l E 



k N 



q=0 



q\ ds<i 



r=2 a=<j 



(4.30) 



s^0 



where 



1 



1 



(2?7 • K) t+1 
(t + l){K 2 Y+ l 



(AT + lJIfofijrfoP+i </ r/)^ 2 ^ (K 2 )^ W U ' 



14 - [* - T^|r,] 



(4.31) 



(r,6,l) 
JV,t 



-1) 



b+1 



P^T^+sfpVT^^ + ai) bKV^) b+1 (Pr,l Pr,2) b ^ 



Pr,l] 



t+1 



t+l(£\K\P r>1 



it+i 



\K + sr, \l) 
|Pr,l - rP r , 2 > 



(4.32) 



6 Notice the summation over r no longer includes r = 1. 
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B 



(r,6,2) 
N.t 



(-1) 



6+1 



d 6 



/j^r^ + s Uy/T=mg. + ai) &Kv^) b+1 (Pr,i p )2 ) 6 *- 6 



Pr, 2 ] 



t+1 



t+l(w r , 2 ] 



t+1 



7V+2 



N+2 



N+2 



+ ST? |£> 

|Pr,2-rP r ,i> 



(4.33) 



5. On pentagon and box formulas 



There are two complexities for box formulas that deserve more study. First, the box formula that we 
discussed before is not a polynomial of u, because it contains also pentagon contributions, indicated by a 
linear factor (ait + b) in the denominator. We need to separate the pentagon part from the box, so that 
the true box coefficient is a polynomial of it. The second complexity is that the null momenta Pj^\{2) (u) 
depend on it in a very nontrivial way (as Qj(u) + x a Qi(u)), unlike the cases of triangles and bubbles. 



For the formulas with a standard input form, these two problems have been solved in 4£, 48 1. We 



review these simplified formulas in Appendix A. 2. Now we want to deal with the general form of tree level 
input. 

In the following two subsections, we first use a "quintuble-cut" method to calculate the pentagon 
coefficients . The true box coefficients can be obtained by subtracting the pentagon contributions. Then 
in the second subsection, we will give a way to simplify the u dependence for the box formula, by generalizing 
the result in the case of standard form. 

5.1 Pentagon coefficient 

The pentagon master integral is 



= J d D P 



• (5-1) 



(p 2 - Af?)((p - K) 2 - M 2 )((p - Ki) 2 - m?)((p - Kj) 2 - m 2 )((p - K r ) 2 - m 2 ) 
The quintuple-cut for the master integral is given by replacing the five propagators with five 5- functions: 

Cut 

d D p 5(p 2 - M't)6((p - Kf - M 2 )5((p - Ki) 2 - mi)5{{p - K 3 f - m 2 )5((p - K r f - m 2 r 



quintuble-cut 

jD m xfJ2 nj2\xff jr , v\2 n J r2\xff„ f.\2 ™2\r//„ f."\2 ™2\r//„ jv- \2 „,2\ 



d 4 £ d~ 2e fi S(£ 2 -n 2 - Mf) 5(-2£ ■ K + K 2 + M 2 — M 2 2 ) S(-2£ -Ki + K 2 + Ml - m 2 ) 



x5{-2£ • Kj + K 2 + Ml - m 2 ) 5{-2£ • K r + K 2 + M\ 



(5.2) 



The integral is totally fixed by five 5-functions. So to get the pentagon coefficient, it is possible to use the 
quintuple-cut method to read the coefficient directly. 



7 The idea of using "quintuble-cut" to determined pentagon coefficients has appeared in 
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As we have discussed before, the coefficient can be accepted as a function of fi 2 . In another word, we 
only need to fix the four-dimensional component of p, i.e. £, and leave the integral for [i. This can be done 
by using the latter four ^-functions. 

The latter four 5-functions give four equations for £ (we have also used the first <5-function that 
I 2 = M 2 + fi 2 ) 



-21 ■ K + K 2 + Ml - Mf 








-2£- Ki + K 2 + M 2 - m 2 
-21- Kj + K 2 + M?-m 2 = 
-21- K r + K 2 + M 2 -m 2 r = 



(5.3) 



by which £ can be solved as 



loK + UKi + LK ) + l r K r , 



(5.4) 



where 



flo\ 

k 

h 
\l r J 



( K 2 



Ki-K Kj ■ K K r ■ K 



Kj -Ki K r - ^ 



K-Ki K 2 
K ■ Kj Ki ■ Kj 



\K -K r Ki- K r Kj ■ K r 



K r ■ Kj 
K 2 



(k 2 



+ M 2 



K 2 + M 2 
K 2 + M 2 
\K 2 + M 2 



M 2 \ 



mi 



m 



m. 



(5.5) 



For a general input T{£) of double-cut integral, the pentagon coefficient can be written as 



Pen[Ki,Kj,K r , K] = T{£) ■ D t {£)Dj{£)D r {£) 



1~(£(i,j, r )) ■ Di(£(i,j,r))Dj(.t(i,j,r)) D r(.t(i,j,r))- 



(5.6) 



In the appendix ^, we show that for the standard tree level input, this formula is equivalent with the 
previous-known formula. We emphasize that the solution £(ij jr ) can a l so De used in generalized unitar- 
ity with multi-cut. For example, for five-cut (with another three cuts across Di,Dj,D r ), we have the 
corresponding pentagon coefficient simply as 



Pen[Ki,Kj,K r ,K] = A\ ree (£) x A^ ee {£) x Af ee (£) x AT e {£) x Af ee (£) 



(5.7) 



There is one point about the solution (|5.4|): there is matrix inverse we need to do. If £(ij iT ) i s contracted 



with any momentum, we may use ( B.4 ) to calculate, which is simpler 
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5.2 Simplify the u dependence of the box formula 



Now we want to simplify the u-dependence of the general box formula. As reviewed in Appendix A..2 , for 
the standard input form, we have known how to simplify the it-dependence of the box formula as 

i / (K 2 r+ 2 UstiiARs 



C[Qi,Qj,K] 



[u 



2 \(£\K\£] n+2 Ut^, (i\Qt(u)\ 

i / (^ r+2 nS(w«w 



W - |Pji,2(«)] + {Pji,l(u) <-» Pji,2(t*)} 

W - l^,i(«)> , 



W - |Pii,a(« = 0)] + {P/i,i(ti = 0) <-► P ii)2 (ti = 0)} 

K) -> |P ii>1 ( u = o)> 

It is worth to emphasize that in the Erst line, u-dependence is everywhere: in the spinor component 
of Pji as well as inside the square roots yf(2Qi{u) x Qj(u)) 2 — AQi(u) 2 qj(u) 2 . Because of this, no matter 
analytically or numerically, it will become very tedious and complicated. To avoid this problem, as analyzed 
in JIgD , another equivalent expression is given as in the second line of above formula, where null vector 
Pji,i(2)(u = 0) as well as square root do not depend on u anymore. The only u-dependence is following 
replacement: 

R s (u) R s (u), Q t (u) Q t (u), (5.8) 

where R s {u) and Qt(u) are given by ( A.20| ) and ( A. 21 ). 

We want to generalize the formula, so that it can be used for the general tree level input. To realize 
this, we need to generalize the above rule, which should not be confined to the special form with R(u) and 



Q(u). This can be achieved if we find a substitution for £. Using ( [A. 20 ) and noticing that ( |AT 

R s {u = 0) 



K 2 



where ot,/3 are given by ( 2.16| ), we can find a relation for R{u) that 



K 2 



{t\K\i] 



\R s \i 



(£\K\£] 



Ps 



(qi,qj,K) 



(a 



(qi,qj,K)s 



( Q («.«)( U ) - l)(-/3g<f - (3P S + (/3 - a)^J^K 



o 



-P 



{£\K\i\ 



[ a (».»)( u ) - l] 



(qi,qj,K) 



■ Pt 



(qi,qj,K)\ 



+ p xx 



1 



K 2 



(/? - a)K 



2£ij ' Ps 



where we have defined 



K 2 



(l\K\t\ 



1] 



(qi,qj,K) 



■ Pt 



(qi,qj,K), 2 



+ P XX 



1 



With the same £ij, it is easy to find that 



K 2 



(£\K\£] 



i\Qt\i 



-2£ij ■ K t + K 2 + Ml 



D 



(5.9) 

I 

■Ps 

(5.10) 
(5.11) 

(5.12) 
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Comparing ( 5,10 ) and ( 5.12| ) with the relation 

K 2 



—2f ■ P 

we can find the rule fl5.8|) is equivalent to 



R s \£], D t {£) 



K 2 



(£\K\£] 



(i\Qt\i] , 



(5.13) 



(5.14) 

(qi,qj,K) 



In Appendix |C], we give another equivalent expression for £ij, which can avoid the appearance of q$ 

Just following the argument of Section 4.2, this new rule can be generalized directly to the formula 
( 4. 14] ) with general tree level input, with only one condition that: £ij, with the substitution for P x j in the 
formulas, should satisfy the two on-shell conditions 



72 



^ 2 - M 2 = 0, (£ij - Kf -fi 2 - Mf = 0. 



(5.15) 



We will proof this is indeed true. 

We first proof the first condition. Using (5.11) directly, we have that (for simplicity we omit the 



superscript of 



ql (£\K\£f &(Qi,Qj) 
where we have used ( |2.14 ) that u = 4/x 2 / (j3 2 K 2 ) , and define 

A(Qi, Qj) = f3 2 {f3 2 [(2 qi ■ q,) 2 - 4q 2 q 2 } + 4^ 2 [a,a j (2 % • Qj ) - a 2 q 2 - a 2 q 2 }} 

By using the relations 8 : 

(P 3hl \K\P Jh2 ] (P ji;2 \K\P ji;l ) = [3 A %[(2 qi ■ qj f - 4q 2 q 2 ] , 

K 2 



(5.16) 



(5.17) 



(P jiA \q \P jit2 ](P ji!2 \K\P jiA ] = 2ip z -pqZJA(Q i ,Q j ) , 



{Pji,i\K\Pji !2 \ (Pji, 2 \qo\Pji,i] 



K 2 

-2^ 2 - 1 -^ A /A(Q ij Q jV . 



we have 



(£\qo\tf 



(£\K\£Y 



VIM 

\e) - l&J) (W\ 2 



'J I 



-4(q 2 ) 2 A(Q u Q, 
l^l^ n(2 qi . qj ) 2 -4q 2 q 2 ] 2 



(5.18) 



Substituting this back to ( |5.16| ), we can find that to have I 2 , — n — M 2 = 0, we only need the relation 



Aq 2 K 2 = (2 qi - qj ) 2 -±q 2 q 2 



(5.19) 



8 These relations can be found in [J49|, below the equation (5.16) in [J49J] . Here we have made a direct generalization to the 
massive case, and also let u — 0. 
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To proof this relation, we notice 



pes 4 



(5.20) 



where p is permutation of fx, 1/2, V3, Using this we can calculate that for q$ 



— w- — 



as 



1 



■1 



(K 2 ) 2 

Noticing that qi ■ k = qj • K = 0, we get 



q 2 K 2 = (q i -q j ) 2 -q 2 ~- 



% Hj ) 



(5.21) 



which is just the relation (|5.19|) . Thus we prove that it's true for the first on-shell condition in ( 5.15| ). 
Then for the second condition, the proof is trivial since we have 

K 2 



2£ii -K = (5- 



%1 " H (i\K\l] x 
Therefore, the general box coefficient is 



(-2P r .K)-(J3- a)K 2 = aK 2 = K l + M{ - M£. 



(5.22) 



C[Qi, Qj, K] 



1 



W -* \Pji,2] + i P ji,l «"* Pjifl] 
\t) - \Pji,l) 



(5.23) 



and the true box coefficient is obtained by subtracting the pentagon contributions: 



Box[Ki,Kj,K} 



)-D. 



E 



Pen.[K it Kj,Kr,lC] 



D 



\e) -» |Pji,i> 



(5.24) 



+ ^ Pji,2} ■ 

Notice that the null vectors -P,j i i(2) are independent of u now. 
5.2.1 Other methods to deal with u-dependence 

In above we have given one way to make u-dependence simpler. In this part we will give other ways to 
deal with n-dependence. 

We can expand any spinor into two independent ones. Sometimes we need to use the spinor or anti- 
spinor components of P\,Pi- For this we can expand into arbitrary null momenta a, b, i.e. Xp = 
and Xp = ^ + p b ^ where t is necessary normalization factor. We can solve that 



(b\a\b] 
(b\P\b\ 



-jam 
(b\p\b\ 



\P) 



-\P\b] (a b) 
y/(b\P\b](b\a\b] 



\P\ 



\P\b)[ab] 



V(b\P\b] (b\a\b] 
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When applying this idea to box coefficients, we use a, b as the Pji^(u = 0) and Pji j2 {u = 0). For the 
box, factor \J (b\P\b\ (b\a\b\ will cancel out eventually as well as factor (a b) , [a b] since numerator and 
denominator are same degree polynomial. Thus we can get final replacement rule as 

\Pji,i(u)) = \P jiA (u)\P jiA (u = 0)], \Pji,i(u)] = \Pji,i(u) \P jitl (u = 0)> , (5.25) 

and 

\Pji,2( u )) = \Pji,2{u)\P jij2 {u = 0)], \Pji,2{u)\ = \Pji,2{u) \Pji, 2 (u = 0)) . (5.26) 



6. Gluon example: A(l ,2+,3 + ,4+,5 + ) 

In this part we use this simple five-gluon example to demonstrate our method 9 . The implementation of 
the formulas into automatic tools is straightforward. In the following example, we do all calculations 



analytically with the Mathematica package S@M [58]. 



By supersymmetric identities [59], the computation is equivalent to one with a scalar field circulating 
in the loop [pCf] . In the sense of four-dimensional unitarity, this amplitude is cut free, i.e. it has only 
rational part and 0(e) contribution. We will use our formulas to calculate whole results including the 0(e) 
contribution, which according to our knowledge, should be the first time. The 0(e) contribution would be 
important in the higher loop calculations, such as discussed in [pi]] , in order to calculate n-parton two-loop 
amplitude to 0(e°), the (n + l)-parton one-loop amplitude needs to be evaluated to 0(e 2 ). 

This five-gluon amplitude can be expanded as a linear combination of one pentagon, five boxes, five 
one-mass triangles, five two-mass triangles and five bubbles. To obtain various coefficients, we need all the 
five kinds of double cuts. We will illustrate the use of the formulas by giving a detail discussion on the 
calculation of the i\23-cut. For the other cuts, we list the final results directly. 

6.1 K 2 3-cut 



The needed tree level input can be obtained from the on-shell recursion relation as [62] 

T K23 (I) = A(£ u 4+,5 + ,l-,£ 2 )A(-e 2 ,2 + ,3 + ,-£ 1 ) 

( (l|^ 2 ^45^i|4] 2 M 2 [4 5] J 



(4 5) (5 1) ((h +fc 4 ) 2 - /i 2 )((£i + K45) 2 -/x 2 )[l|i^l|4] Kf 51 [5 1] [l|ir 45 *i |4] 
M 2 [2 3] 



(2 2,){{h-k z f-^) t 
H 2 [2 3] (l|(if 23 - £)K A5 £\A} 2 //[2 3] [4 5] ; 



(4 5) (5 1) (2 3) D 1 {£)D 2 (£)D 3 {£) [1\K 45 £\A] ^ 45 i ( 2 3) [5 1] Dt^K^] 



(6.1) 



9 We emphasize that the main purpose of choosing this simple example is to illustrate the using of our formulas. Since the 
tree level inputs will have spurious poles, it can also serve as a good example to check our formulas. The same procedure 
should be implemented directly to compute more complicated cases. 
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where we have defined 

£ = h = K 23 -£ 2 , (6.2) 

and 

K = K 23 , K 1 = k 3 , K 2 = -k 4 , K 3 = -K 45 , (6.3) 

Di(i) = (3\£\3], D 2 (£) = -(Mm, D 3 (I)=Kl 5 + 2l-K 45 . (6.4) 
Notice that there is a spurious pole in the tree level input: 

[11^45^4] = [1 4](4|F|4] + [1 5](5|£j4], (6.5) 



so it can serve as good example for checking our generalized formulas. The first term of the input ( |6.1[ ) 
has degree 0, while the second term has degree —2, which does not contribute to triangle and bubble. 
Furthermore, by comparing ( |6.4[) and (|6.5]), we can find that there's no subtle relation between the spurious 



pole and D 2 {£) or D 3 (£), as that mentioned in Section 3.4, so this second term has no box (and no pentagon) 



contributions either. Therefore, we only need to consider the first term in ( |6.lD in the calculation, i.e. we 
can let 

T (J) - ^ 3] (11(^23-^45^]^ 

Km{ > (4 5) (5 1) (2 3) Dl {i)D 2 {£)D 3 {£) [1\K 45 I\A] ' 1 ' ' 

From this cut we can calculate the coefficients of the pentagon, three boxes, one one-mass triangle, 
two two-mass triangles, and one bubble, as shown in Figure 2. We calculate them case by case. All the 
needed formulas are collected in Section [T~l| . 

6.1.1 Coefficients of pentagon 

There is only one pentagon. Using ( |i~TT|) , we have 

Ven[K ly K 2 ,K 3 ,K] = \r K23 (t) ■ D X {1)D 2 (£)D 3 (£) 

^^(1,2,3) 

^ 2 [2 3] (l|(^23-Al,2,3))^45Al,2,3)|4] 2 



(4 5) (5 1) (2 3) [11^(1,2,3)14] 
s 3 23 sl 5 s 12 s 15 s u (1 2) (3 4) (1 5) 2 



A 3 

where A is defined as 



M (6-7) 



A = (2 3) ((4|fcifc2*afci|5> + {4\kihk 2 k 3 \5) + (4|fc 2 fe 3 A; 4 A;i|5>) . (6.8) 
Obviously, the pentagon contribution is 0(e). 
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1/ l-K 3 \5 
Pen[K 1 ,K 2 ,K s ,K] 



Box[K!,K 2 ,K] 



2 / 1 M 
Boxfi^i, K 3 , K] 



Box[K 2 ,K 3 ,K] 




Tri[K u K] 






Tri[K 2 ,K] 



Tri[K 3 ,K] 



Bub [A'] 



Figure 2: The first figure illustrates the ^23-cut. And the other figures show all the scalar bases of which the 
coefficients can be obtained from this cut. 

6.1.2 Coefficients of box 

There are three boxes. Using fll.l2p , the box from K\,K2 is 
c [5i|2|3|4] = Box[Ki,K 2 ,K] 



r K33 (e 12 ) ■ D 1 (e 12 )D 2 {e 



12 



Pentagon 



\P\ -* |P2l, 2 ] + {^21,1 <~> 1*21,2} 
W - |ftl,l> 



D 3 {£ 12 ) _ 

(4 5)((1|2|4|3|2|1|5]+(1|3|2|4|3|1|5]+(1|5|4|2|3|4|5]) 2 . 2 . 2 , , , 



si 5 s 23 2 (2 4) (3 4) A 
where we have changed u to jti 2 , by using the relation 



-(H z y + {n -term) 



u 



4fi 2 
K 2 ' 

^ 23 



Similarly, we can get another two box coefficients 



C[23 |4|5|i] = Box[K 2 ,^ 3 ,K] = MMMM! (/i 2 )2 + (/i 2. tenil) 



A 



c [45|1|2|3] = Box^,^] = (1 2)(1 3 ^ 2]2[5 V ) 2 + (/Aterm) 



(6.9) 



(6.10) 



.11) 



(6.12) 



As expected, after subtracting the pentagon contribution, the box coefficient is a polynomial of fi 2 . The 
complete box coefficients with ^ 2 -term will be given in Appendix 
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6.1.3 Coefficients of triangle 

There are three triangles. Using with N = 0, the triangle of K\ is 



3 [451|2|3] = Tri^i,^] 



2 (-y/T=u) yJ—Aq{K 2 (Pi,i Pi,2 
d ( {£\K\t\ 



dr \ (K 2 ) 



W - \Qi \t) + {Pi,i 

\l) \Pi,i - rPi, 2 > 



((1 4) 2 (2 3) 2 -2(l 3) 2 (2 4) 2 )(1|2|3] 2 
2(1 5) (2 3) 2 (2 4) 2 (3 4) (4 5) M 

Similarly, we get another two triangles 

— Trl\K K\- (1 2)(1 4) 2 (4|K 23 |4] 
c [4 | 51 | 23] - m[K 2 ,K\ - " 2(1 5)(2 3)(2 4)2(3 4)(4 5) 

C[i[23|45] = Tri[K 3 ,i^] = 

6.1.4 Coefficient of bubble 

Using ( |1.14| ) for n = and A; = 3, we have 



Bub [if] = if 2 



r=l 



M,l), 



where 



T(£) 



K> - \K\ V ] 



B 



(r,0,l) 
0,0 



(f,0,2) 
0,0 



(0) 



(0) 



2r/-iT 1 
K 2 (Z^) 

i i f (e\K\£] 

VT^l y/=4%K* (£\K\P r)1 ] (e\ V K\£) V (i^ 2 ) 

1 <%|Pr- )2 ] 1 f(AK\£] 



u 



^Z4gK2(e\K\p ri2 ](e\r,K\e) V (k 2 ) 

The explicit results for each terms are 

(l\K 23 |r?] 3 



• D r (£) 



<S(0) 



B 



(1,0,1) 
0,0 



B 



(1,0,2) 
0,0 



B 



(2,0,1) 
0,0 



(0) 



(0) 



(0) 



S23<1 5) (2 3) 3 (4 5) [2 r/][3 77] <4| i^ 23 
(1 3) 3 [3 t?] 



s 23 (l 5) (2 3)3(3 4) (4 5) [2 rj\ 

(1 2) 3 [2 77] 

S23 <1 5) (2 3)3(2 4)(4 5>[3 77] 

(1 4) 3 (5|4|7y] 

(1 5) (2 3) (2 4) (3 4) (4 5) 2 (4|if 23 |4](4| J FC 23 |7 ? ] 



r^0 



(6.13) 

(6.14) 
(6.15) 



(6.16) 



w - 




w - 


- \Pr,l) 


14 - 




l«> - 


- |Pr,a) 



(6.17) 
(6.18) 
(6.19) 
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and 

£gf 2) (0) =^ 3 °' 1) (0)=^ 3 °' 2) (0)=0. 

After taking them back in ( 6.16[) and doing the summation, we can find 

(1 2) 3 (3 4) 2 [4 2] + (1 3) 3 (2 4) 2 [4 3] - 3(1 2)(1 3)(1 5)(2 4)(3 4) [5 4] 2 
c [23 |45i] = Bub [A 23 ] = (1 5)(2 3)3(2 4}(3 4)(4 5)(4|A 23 |4] (6 - 20) 

As expected, r/ doesn't appear in the final result. 
6.2 Other cuts 

For other cuts, we first give the necessary tree level input for the calculation |3|], then list the final results 
directly. 

6.2.1 Tree level input 

The needed tree level input obtained from the on-shell recursion relation for K^-cut is 

T Ks4 (I) = A(e 1 ,5 + ,l~X,h)A(-e 2 ,3 + ,4+,-£ 1 ) 

( (1|^|5] 2 (1|^|2] 2 ^ 2 [5 2] 4 \ 



(5 2)((£ 1 + fe5) 2 - / u 2 )(^i + ^5i) 2 -/x 2 )[2|K5i^i|5] K| 12 [5 1][1 2][2|K 51 4|5] 
^ 2 [3 4] 



(3 A)((£ 1 -k 4 y-^) / 
^ 2 [3 4] (l\7\5]i(l\(K u -I)\2] 



(3 4) (5 1) (1 2) ((£ _ kA )2 _ ^ W + h) 2 _ M 2) ((i + k 51 )2 _ M 2) [ 2 |iT 51 £|5] 

_ //[3 4] [5 2] 4 1 

~ #512 ( 3 4 ) I 5 !][! 2] ' ((£_ fc 4 )2 - ^)[2\K 51 I\5] 
and -ftT 2 3-cut it is 

T Kl2 (I) = A(£ u 3 + ,A + ,5 + ,e 2 )A(-e 2 ,l~,2 + ,-h) 



(6.21) 



(3 4)(4 5)((£x + A;3) 2 -^ 2 )(^i+^34) 2 -^ 2 )y ^ K\ 2 ^l x - k 2 f - ^) ^ 
M 2 [5\K U 7\3] (1|£|2] 2 



(3 4) (4 5) K\ 2 ((£ _ fc 2 )2 _ /U 2 )(( ^ + ^2 _ M 2)((£ + K34 )2 _ ^2) ' 



(6.22) 



Similar to the case of if 2 3-cut, the second term in Tk 34 (£) can also be neglected. ^45-cut and ifsi-cut are 
not necessary, since they are related to X 2 3-cut and i^i 2 -cut by symmetry, although we have also done 
the calculation independently in order to show the results are consistent. The coefficients are related by a 
symmetry operation, for example 



c [51|2|3|4] = -C[12|3|4|5] 
which is similar for other coefficients. 



(6.23) 

2^5,3^4 
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6.2.2 Other coefficients 

The other two boxes are 

(2 3)((1|2|3|5|4|3|2] + <1|4|5|3|4|1|2] + (1|5|3|4|5|1|2]) 2 2 2 2 
C[ i2|3|4|5] = , 12S45 2 (3 4)(3 5)A ^ ] + ^ " term) (6 - 24) 

_ (1 2)(1 5) (2 3) (4 5) [4 3] [5 2] 2 2 2 , 
C[34|5|i|2] - { 2 5^3 4 ^ A (/* ) + (/* " term ) ( 6 - 25 ) 

The other four one-mass triangles are 

(1 2) 2 (1 5)[2|1] , 

^l 1 ! 2 ! " 2(2 3)(2 5) 2 (3 4)(4 hf (6 - 26) 

_ (2(1 4) 2 (3 5) 2 -(l 3) 2 (4 5) 2 )(1|5|4] 2 
C [ 123 I 4 I 5 1 " 2(1 2) (2 3) (3 4) (3 5) 2 (4 5) 2 ^ (6 - 27) 

- < 1 2)(1 5)2[5 1] , 2 rfi 28^ 

^l 5 ! 1 ! " 2(2 3)(2 5) 2 (3 4)(4 hf (6 - 28) 

= [f_3j , . 

C[5i2|3|4] 2(1 2)(1 5)(2 3)(2 4) 2 (3 4) 2 (3 5) 2 (4 5) 3 1 j 

((1 4) 4 (3 5) 2 ((2 5) 2 (3 4) 2 - 2(2 4) 2 (3 5) 2 ) - 2(1 4)(1 5) 3 (2 4) 2 (3 4) 3 (3 5) 
+2(1 4) 3 (1 5)(2 4)(3 4)(3 5) 2 ((2 3) (4 5) + (2 4)(3 5)) + (1 5) 4 (2 4) 2 (3 4) 4 )/x 2 

The other three two-mass triangles are 

_ (1 5) 2 ((1 3)(2 5) + (l 2)(3 5))<5|g 12 |5] 2 
C 1 5 I 12 I 34 ] " 2(1 2)(2 5) 2 (3 4)(3 5) 2 (4 5) (6 - 30) 

_ (1 2) 2 ((1 5)(2 4) + (l 4)(2 5»<2|Jr 51 |2] 2 
C 1 2 I 34 I 51 ] " 2(1 5) (2 3) (2 4) 2 (2 5) 2 (3 4) (6 - 31) 

(1 3) 2 (1 5)((3|4|3] + (3|5|3]) 2 
C 1 3 I 45 I 12 ] " " 2(1 2) (2 3) (3 4) (3 5) 2 (4 5) (6 - 32) 

The other four bubbles are 

-3(1 2)(1 4)(1 5)(3 4)(3 5)[3 2] + (1 4)3(3 5) 2 [4 3] + (1 5) 3 (3 4) 2 [5 3] 2 
C[45|123] " (1 2)(2 3)(3 4)(3 5)(4 5)3(3|K 45 |3] ^ [ ' 

f (l 3) 3 ((2|1|5]((3 5)((2|4|1] + 2(2|5|1]) + (2 5)(3|4|1]) - (2 3)(2|4|5|2|1]) 
[12|345] " I (2 3) 3 (3 4)(3 5)(4 5)«2|1|5|3) - (2|5|4|3))(3|K 45 13] 1 " ' 

(1 2) (£i5(l 3) (2 3) + 2(1 3)(2|5|2|3) - (1 2)(3|4|5|3)) \ 2 
(2 3) 2 (2 5)(3 4)(4 5)((2|1|5|3)-(2|5|4|3)) / 

/ (l 4) 3 ((5[1|2]((2 5)(4|3|1] + (2 4)(2(5|2|1] + (5|3|1])) - (4 5)(5|3|2|5|1]) 
C[51|234] " V (2 3)(2 4)(3|4)(4 5)3(-(4|2|1|5) + (4|3|2|5))(4|K 23 |4] l " > 

(1 5) (-512(1 4) (4 5) + (1 5)(4|2|3|4) - 2(1 4)(4|5|2|5)) \ 2 
(2 3)(2 5)(3 4)(4 5) 2 (-(4|2|l|5) + (4|3|2|5)) ^ 



- 35 - 



(1 5) 2 (3 4)(1|5|1|5)(3|4|3|4) 



| (1 3) 2 (1|3|1|3) | (1 4) 2 (1|4|1|4) 



c [3 4|5i2] (1 2)(1 5)(3 4)3 5)(3 5)(4 5)( 5 |2|i](5|^ 12 |5] 1 (2 3) (3 5) [4 1] ' (2 4) (4 5) [3 1] 
((4|1|2|5) - (4|3|1|5)) 4 , (1 2) 4 (3 4)»«5|1|2][3 1] - <5|2|1][3 2]) 



+ 



(4 5)3(2|5|1](5|2|1](5|4|1][3 1] (2 3)(2 4)(2 5)(-(4|2|l|5) + (4|3|2|5))(2|5|1] 



As it is obvious, one given coefficient can be calculated from various cuts. Using Mathematica we 
calculate them and do find same results. Also, we check that the coefficients satisfy the symmetrical 
relation ( 6.23| ). These are strong tests for our method. 

6.3 Rational parts 

To extract the rational parts, we need the following relation [50, 42 1 



+ 0{e), /IV 1 



+ 0(e), jf[// 2 ] 



K 2 



+ 0(e) 



The rational part of the amplitude is 

1 



R 



(47T) 



6(/i 



2\2 



c [12|3|4|5] + c [23|4|5|l] + c [34|5|l|2] + c [45|l|2|3] + c [51|2|3|4]) 



(/i 2 ) 2 -term 



1 



+ ^2 l C [123|4|5] + c [234|5|l] + c [345|l|2] + c [451|2|3] + c [512|3|4] 



+ c [12|34|5] + c [23|45|l] + c [34|51|2] + C[45|12|3] + c [51|23|4], 



1 



2 2 2 2 2 \ 

- 6^2 ( C [12|345]-^12 + c [23|451]-^23 + c [34|512] -^34 + c [45|123] K 45 + c [51|234]^51y 



The rational part of the same amplitude has been given in 
i 1 



as 



R 



48vr 2 [1 2] (2 3) (3 4) (4 5) [5 1] 



S23 + s 3 4 + s 45 )[2 5]" - [2 4] (4 3) [3 5] [2 5] 



1 2] [1 5] ( 2 2 [23] + [34] + 2 2 [45] 



(12) (15) V N ' (2 3) * ' * ' (3 4) 

With the help of Mathematica one can easily check the following result 

r=\r 



(6.37) 



(6.38) 



(6.39) 



(6.40) 



The factor 1/2 is because R is for the amplitude with a real scalar circulating in the loop, while that for R 
it is a complex scalar. As we can see, the spurious poles, such as A in the box coefficient, are all cancelled. 

Our result R looks more complicated than the R. The reason is because when we use the Mathematica, 
we have chosen a basis for spinors and expressed other quantities by this basis. The choice will sacrifice 
some simplicity, but make the result unique for the chosen basis and easy to compare. Furthermore, it 
makes the evaluation straightforward either analytically or numerically. 
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7. Summary 



In this paper, we generalize formulas for coefficients obtained by D-dimensional unitarity method Eg, 47, 



48, 49 1 to the case where spurious poles are allowed in tree-level input. While keeping algebraic expressions 
explicitly, this generalization makes sure that we can use the most compact tree-level input to reduce the 
complexity of computations. We summarize some main points of our final formulas as following: 

• The tree-level input (obtained by any method) can be inserted into the formulas directly without 
any modification. 

• The formulas are for the full coefficients, i.e. including all the contribution of cut part, rational part 
and 0(e) part. If we want only the cut part, we can take the four-dimensional tree level input, and 
let u = in the formulas. This will simplify the formulas dramatically. 

• The formulas are for the general case, massive or massless. 

• The formulas are at the analytical level. It can give compact analytical expression for some results, 
such as that shown in the example of Section 6. 

• The formulas are suitable for the non-renormalizable case, such as the one loop calculation of gravity 
theory. 
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A. Review of formulas 

We review formulas of various coefficients for the standard input form |46|, ^7|, |4^, ^| 

r (n),a_ n?^-2i-p 3 o _ cgr ngfwij (AU 
sf[) nix aw (wrnUmw [ ' } 

We have defined the quantities: 

Rj{u) = -p{VT^l) Pj + pjK, Qi(u) = -p(y/T=u) qi + mK, (A.2) 

where 

Pa ■ K 



Pi = [PJ- 

K, • A' 



= K, 



^_<^( 1 + M_M), (A ,) 

4 «=-^(l + ^) + £±^l (A,) 
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and (5 are given by (|2.16|) . 



A.l Formulas 

The box coefficients (pentagons not separated) are given as 

{K 2 )n+ 2 f na^m 



C[Qi,Q v K] 



2 \(p^i\m*M n+2 nli, t ^ {Pji,i\Qt\Pjid 



+ {P jitl «- P jiy2 } , (A.5) 



where 



2^2 



Aji = (2Qj ■ QiY - \Q]Q 
Pji,i = Qj + ( — ® J ^ 9 ) Q, 



Pji,2 — Qj + 

The triangle coefficients are given as 
{K 2 ) n+1 1 



-2Qj • Qi -y/Ajj \ 

2« r j 



C[Q S ,K] 



n+l 



2 ( v / Al)" +1 (n + l)!(P s , 1 P Si2 

/ (P^-rPs^RjQslPs^-TPs^) 



■ dT n+l \j[* =1 ^ s (P Sjl - TP s ,2\QtQs\Ps,l ~ tP s , 2 ) 



T = 



where 



(A.6) 



(A.7) 



A s = (2Q S • A-) 2 - 4Q 2 A 2 

Note that the triangle coefficient is present only when n > — 1. 
The bubble coefficients are given as 



C[K] = {K 2 ) n+1 ^ 



n+1 » (_!)„ d9 / 



9=0 



r=l a=g 



C-^)+EEK- ;1) ( s )-<- 2, ( s 



,(r;a-g;2) , 



where 



_ 1 (2r / -A)*+ 1 ngfm(^+^7Hg 

dr" y n \[ri\rjK\rj\» (t + 1)(A 2 )*+! (I r,) n+1 rg =1 (£\Q p (K + s V )\£) 



\\t)-+\K-Trj\ V ] 



(A.8) 



(A.9) 



(A.10) 



T = 
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(-1) 



6+1 



1 (Pr,l ~ TPrMv\Pr,l] t+1 

b\^ h+1 (P rsl P r , 2 ) b dr b {(t + l) (P r>1 - TP r , 2 |A|P r>1 ]' +1 

rn+fc 



(P rjl - TP r , 2 \Q rV \P rA - rP r , 2 >° n"=r <-Pr,i - rP r ,2\Rj(K + srj)\P rA - T P, 



r,2 



(P r ,i - rP r , 2 \ V K\P rA - rP r , 2 ) n+1 Up=i^ r ( p r,i ~ rP r , 2 \Q P (K + s V )\P rA - rP r , 2 > 



(A.11) 



T=0 



„(r;b;2), , d»( 1 <P r , 2 - rP rA MP,, 2 ] m 



blVK b+1 (P r>1 P r , 2 ) b dr b \(t + l) (P r , 2 -TP r ^\K\P^ 2 } t+1 

(i>r,2-TP r ,l|Q r 77|P r , 2 -TP r ,^ \ 

x ± ] { P± 12) 

(Pr,2 - rP^ilr/^lP^ - rP ril )™ +1 nJU^r <^,2 - rP r ,i|Q p (A + ^)|P P)2 - rP r ,i) / t=q ' 

where A r , P r) i, P r>2 are given by ( |A.8| ), and 77,77 are arbitrary, generically chosen null vectors. Note that 
the bubble coefficient exists only when n > 0. 

A. 2 Formulas with u simplified 

The u-dependence for above formulas can be simplified further ||49| , |48| . Here we collect these simplified 
formulas. The pentagon formula is also given. The most important simplification is for the null spinors. As 
we have emphasized in main text, in the simplified formulas, u-dependence becomes minimum and much 
simpler, especially the null momenta do not depend on u anymore.. The null spinors Pji t a (o = 1>2) in 
box formula are taken as 

\Pji,a) = \Pji,ai u = °)> . \ P iiA = \Pji,a( u = °)l ( A - 13 ) 
while the null spinors P SiQ in triangle and bubble formulas are taken as, 

\Ps,a) = \P qs ,a) , \Ps,a] = \P qs ,a], (A.14) 

where 



P qa ,a = q s ± I \K. (A.15) 

Notice that for simplicity we use the same symbols for the null spinors as those in Appendix |A1| Whether 
they depend on u or not should be clear according to the context. We also use different notations for the 
coefficients, such as Box[Qi,Qj,K],Tri[Q s ,K] other than C[Qi, Qj, K], C[Q S , K]. 
The pentagon coefficients are given by 

Pen[Q uQj ,Q t ,K] = (A 2 )" +3 lla=1 * (A.16) 
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where 

J Ki ,K 3 .,K t ;P B ) _ (Kf + M 2 - m 2 )e(P s , Kj, K, K t ) + (K 2 + M 2 - m 2 )e(K t , P s , K, K t 



K 2 e(K i} Kj,K, K t ) 

{K 2 + M 2 _ M 2 )e{K h Kj, P s , K t ) + (K 2 + M 2 - m 2 )e(K t , Kj,K, P s ) {U _ 



K 2 e{Ki, Kj,K, K t ) 



(K^M) = (Kj + Ml - m 2 )e(K, Kj,K s , K t ) + (K 2 + M 2 - m 2 )e(i^ K, K s ,K t ) 
ls K\(K t ,Kj,K,K^ 



(K 2 + M 2 - m*)e(Ki, K 3 ,K, K t ) + (K 2 + M 2 - m 2 )e{Ki, Kj,K s , K) 



K 2 e(Ki, Kj ,K,K t ) 
{K 2 + M 2 - M 2 )e{K u K 3 ,K S , K t ) 



K 2 e(Ki, Kj,K, K t ) 
The box coefficients are given by 



Bax.[Qi,Qj,K\ 



2\2+n 



l& (P j i i l\R.(u)\P jifi 



{ (P^i\K\p Jt;2 ] n+2 u k t=1 , t ^ (Pji-MQtiu)]^ 



E 



t=l,t^i,j nt=l,t«^i,j,t 7w " " "" ^ {Pji;l\Qt(u)\Pji;2 



(A.18) 



(A.19) 



where 



and 



(qi,qj,K) 



Rs(u) = - ( a («^) (tt) _ + jRs(u = o), 



lqi,qj,K) 



(9o * ) 



K 2 



(A.20) 
(A.21) 



(A.22) 



(u) = 



/3 2 (l-u) + 



4K 2 [a i a j (2q i -q j )-a?q?l-a 2 j q?] 
(2gi-gj) 2 -4g^g^ 



/3 2 + 



4JiT2[a i a J -(2 qi - gj -)-a| g |-a^ t ?] 
(2gi-</i) 2 -4</?d 



(A.23) 
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The triangle coefficients are given by 
{K 2 ) 1+n 1 



TrifQs, K] 



2 (-/VT=^rH-i( v '=4^K2)»+i (n + 1)! <P S)1 P s , 2 ) n+1 
(T+i / nS (Ps,i ~ TP e , 2 \Rj(u)Q a (u)\P ail - tP 



s,2 



dT n+l (Ps)1 _ rP s>2 \Q t {u)Q s {u)\P s ^ - rP s , 2 ) 

The bubble coefficients are given by 



+ {Ps,l «- Ps,2} 



(A.24) 



T = 



1 d<? 



Bub[K] = (K 2 ) i+n y: ~ I agUw «- a 9;1) ( S ) - <rr 2) « 



g=0 



(A.25) 



s=0 



where 



dr« I nlfal^M" (t + (£ ^"+1 fj* (£\Q p ( u )(K - srj)\i) 



,(A.26) 



r=0 



(-i) 



6+1 



1 (P r ,i - rP r>2 |T;|P r) i] 



t+i 



6!( _ / VT^)myZ4^2 b+1 (p p ^dr*- ^(t + 1) (P rA - rP r , 2 \K\P rA ] t+1 



(P r ,i - TP r , 2 \Q r (u)ri\P r ,i - rP r , 2 ) b nS (^r,i - TP rt2 \Rj(u)(K - s V )\P rA - tP, 



r,2 



(P r ,x - rP r , 2 \r,K\P rA - rP r>2 ) n+1 U P =i, P ^ r ( p r,l ~ TP r , 2 \Q P (u)(K - srj)\P r ,i - rP r , 2 ) 



,(A.27) 



T = 



■1) 



6+1 



t+1 



1 {P r ,2 ~ TP r M P r,2 

t+1 



b\(-pVT^) b +W-^K 2b+1 (P r>1 P, 2 ) b ^ + (P;2 - rP rjl \K\P r , 2 



(P r , 2 - TP r ,i\Q r (u)ri\P r , 2 - rP r ,i) 6 nS ( p r2 ~ T P r,l\Rj{u){K - sr,)\P^ 2 - TP r ,l) 



(P r ,2 - rP rA \ V K\P rj2 - rP ril ) n+1 Up=i, P ^ r ( p r,2 ~ rP r ^\Q p (u)(K - s V )\P r , 2 - rP, 



.(A.28) 



T=0 



All the u-dependence is explicitly presented in the formulas. They are only from R(u), Q(u), R(u), Q(u) 
and the yl — u factors. 



B. Compare the pentagon formulas 

For the standard form input 

T(7) 



n:i fe (-2^-p) 



the pentagon coefficient by using (|5.6[) is 



Utl(-^(iJ,r)- P s 



Pen[Ki,Kj,K r ,K] 



ULx, m)j ,r(- 2i (i,j,r) ■ K * + K t + M l " m t) 



(B.l) 



(B.2) 
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where %j,r) is given by (5.4). We need to show that the above formula is equivalent to the pentagon 

jn+k a {Ki,K h K r -P s ) 



coefficient flA.lq ) 



-2\n+3 lls=l Pb 



Pen[Qi,Qj,Q r ,K} = (K z ) n ^ {Ki , Kj , Kt , Kr y 



(B.3) 



where p[ K " K i> K « p °) ^ ^Kije jt K tt K r ) are giyen by and flPg|) . 

To prove the equivalence, we first show that 



2 o(Qi,<lj,1r;ps) 



By using the expansion for P s as 



where 



2£ {iij>r) ■ P s = K z p, 



P s = b K + biKi + bjKj + b r K r 



(B.4) 
(B.5) 



fb \ ( K 2 Ki ■ K Kj ■ K K r ■ k\ ( K ■ P„\ 



.1 

\b r J 



KKi K 2 Kj -Ki K r - Ki 
K ■ Kj Ki ■ Kj Kj K r ■ Kj 
\K ■ K r Ki- K r Kj ■ K r Kj J 



Ki-P s 
K r P s 
\K r ■ P s ) 



e(K h Kj,K,K r 



fe(K h Kj,P s ,K r )\ 
e(P s ,Kj,K,K r : 
e(Ki,P 8 ,K,K r ) 

\e(Ki,Kj,K,P s ) J 



and Q5.5| ), we can find that 

-2(l K ■ P s + hKi ■ P s + IjKj ■ P s + l r K r ■ P s ) = (K 2 + Ml - M 2 2 )6 + (K 2 + M\ - m 2 )bi 

+ (K 2 + Ml - m 2 )bj + (K 2 + Mf - m 2 r )b r . (B.6) 



This is just the relation (B.4). Then, with this relation, we can find directly that 

-2? (Wfr) • K t + K 2 + Ml -m 2 = K ^K h K t ,K r) (B ?) 

With (|B.4|) and (|B.7|), it's obvious that the two pentagon formulas ([B.2D and QB.3]) are equivalent. 



(qi,qj,K) 



C. Another equivalent expression for 



As in section we may give another equivalent expression for £{j that avoid the appearance of q 
By using the expansion 



Ps = a 
where the coefficients are: 



(qi,qj,K;p s ) (qi,qj,K) (qi,qj,K;p s ) , (qi,qj,K;p s ) 



+ a 



gi + a) 



gj> 



a, 



_ {Ps-qt ,qi,K) ) _ e(P s ,K u K 3 ,K) 







(quqj,K)\ 



K 2 (Q K o 



(qi,qj,K), 



{qi,q h K; Ps ) _ (Ps ' gi)gj ~ (Pa ; gj)(gi ; gj) _ J quqj ) J qi , qj ) _ g*gj ~ gj (gi " gj) 



g? - (gi • gj) S 



{ qi ,qj,K;p s ) _ (Ps • gj)g 2 ~ (Ps ■ gi)(gi • gj) 



^ 2 g? - (gi • gj) 5 



gfgf ~ (n ■ gi) 2 ' 

Ps . tfe*) = ml - n(n ■ qj) 

3 gfgj - (gi • gj) 2 ' 



(C.l) 

(C.2) 

(C.3) 
(C.4) 
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we can write 



R s (u) 



Ps ■ q 



(gi,gj,K)s 



- l)(-/3g(*'*'^) + i? s ( u = 0) 



_ i) (a\ 9i ' qj ' K;Ps) Qi(u = 0) + = 0) + A 



(qi,qj,K;p s ) 



K - R s (u = 0) 



+ P s (u = 0) 



(gi,gj,K ;p s ) 

-(a^\u) - i)p^'is> K '^K + a^*)(u)i? s (n = 0), 
where we have defined 



(C.5) 



Q: ; 



(g%,gj) 



P 



(C.6) 



(P^i(u = 0)|Q,(u = 0)|P iW (« = 0)] = 0, 



Making use of the properties 

{P jiA (u = 0)\Qi(u = 0)\Pji, 2 (u = 0) 
we have 

R s ( u ) = -( a («-«)(«) - i^^^k + a^\u)R s (u = 0). 
Using (p|) and flClj ), we can further rewrite it as 



(C.7) 



R s (u) = (a^\u) - 1) 



K 



(9i>9j) 



P., 



K + a(**J(«) f-/3P s + (/3 - a) 



P,-K 



K 



(a<«'«>(u)/3 - a)^(K ■ P s ) + (a(*'^(«) - 1) [(a^ fe ' %) + 



Ps 



K 2 

K -a^\u)pP, 



Then we have 
K 2 



(£\K\£] 



t\R s \£ 



l r 

K 2 



-a 



(<2i 



Therefore we find the equivalent rule is that 



^(u)p XJ y Ps 



1 

2 



<l>) 



{£\K\£j 



(C.8) 
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D. Complete Box coefficients of the example 

In this appendix, we give complete box coefficients. 



C[23|4|5|l] = 




(D.l) 



C[45|l|2|3] - 




(D.2) 



<4|5>«1|2|4|3|2|1|5] + <1|3|2|4|3|1|5] + (1|5|4|2|3|4|5]) 2 , 2 , 2 , (1|2)(2|3) 3 [3|2][4|3] 



{3(1|4) 2 (1|5)(4|5) 2 (1|4|1|4)(4|2|3|4)(-(5|1|3|2|1]-(5|3|2|5|1]) 2 

+ (1|4) 5 (4|5) 2 (-(5|1|3|2|1] - (5|3|2|5|1]) 3 - 2(1|5)(4|5)(1|4|1|4)(1|4|1|5)(4|2|3|4) 2 

«5|1|2]<4|5)(-<2|3|1] + <2|5|1]) + 2(4|5)(5|3|2|5|1]) 

-2(1|5) 2 (1|4|1|4)(4|2|3|4) 2 ((4|5)(-(2|5)(4|3|1] + (2|4)(-3(5|2|1] - 2<5|3|1]))<5|1|5|1|2] 
-(4|5) 2 (5|3|2|3|2|5|1] - (1|5) (2(2|4)(3|5)(4|2|1](2(5|2|1] + <5|3|1]) 

-(2|5)(3|4)(4|3|1]((5|2|1] - <5|3|1]) + (2|4)(3|4)(-2(5|2|1] 2 + 3(5|2|1](5|3|1] + < 5| 3 1 1] 2 ) ) [3| 2] ) 
+ (1|5) 3 (2|4)(4|5)(4|2|3|4) 2 (4|^ 51 |4] 2 (4|5)([2|1|5|1] - [2|3|2|1] - [2|4|5|1]) 
+ <1|4><1|5> 2 <4|2|3|4> 2 (<4|5>(<2|5><4|3|1] + <2|4> (6<5| 2| 1] + 5 <5 1 3 1 1] ) ) < 5| 1 1 5| 1 1 5 1 1 1 2] 
-<4|5> 2 <5|3|2|3|2|3|2|5|1] + (1|5|1|5)( - (2|4)(3|5)(4|2|1](-13(5|2|1] - 10(5|3|1]) 
-(2|5)(3|4)(4|3|1](2(5|2|1] - (5|3|1]) - 2(2|4)(3|4)(5(5|2|1] 2 + (5|2|1](5|3|1] - <5|3|1] 2 )) [3|2] 
-(1|5)( - (2|4)(3|5)(4|2|1](-8(5|2|1] - 5(5|3|1]) - (2|5)(3|4)(4|3|1]((5|2|1] - 2(5|3|1]) 
+2(2|4)(3|4)(-2(5|2|1] 2 + 2(5|2|1](5|3|1]+(5|3|1] 2 ))[3|2|3|2])}^ 2 (D.3) 



C[51|2|3|4] 



2 (2|4)(3|4)A 



X 



(1|5)(2|4)2(4|5)3 
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_ <2|3)((1|2|3|5|4|3|2] + (1|4|5|3|4|1|2] + (1|5|3|4|5|1|2]) 2 2 2 _ (1|5)(4|5) 3 [4|3] [5|4] 
[12|3|4|5] " s 12 s 4 5 2 (3|4)(3|5)A ^ ' 2A 3 (1|2)(2|3) 3 (3|5) 2 

{2(1|3)(2|3)(1|3|1|2) 2 (3|4|5|3) 2 (2(2|4|5|2)(3|2|1] + (1|3) 5 (2|3) 2 ((2|1|5|4|1] + <2|5|4|2|1]) 3 

-3(1|3) 3 (2|3) 2 (1|3|1|2)(3|4|5|3)((2|1|5|4|1] + (2|5|4|2|1]) 2 + (2|1|5](2|3)((5|4|1] - <5|2|1]) 

-2(1|2)(1|3)(1|3|1|2)(3|4|5|3) 2 ((2|3)(2|4|5|4|5|2)(3|2|1] - (2|3)((3|5)(2(2|4|1] + 3(2|5|1]) 

+<2|5><3|4|1])<2|1|2|1|5] + < 1 1 2> ( <3 1 4> <3 1 5> ( <2 1 4 1 1] 2 -K 3 <2| 4| 1] <2 1 5 1 1] - 2<2|5|1] 2 ) 

-<2|5)<3|4)(-<2|4|1] +(2|5|1])(3|4|1] - 2(2|4)(3|5)(-(2|4|1] - 2(2|5|1])(3|5|1]) [5|4]) 

+(1|2) 3 (2|3)(3|5)(3|4|5|3) 2 (3|K 12 |3] 2 (2|3)([5|4|5|1] + [5|3|2|1] - [5|1|2|1]) 

+(1|2) 2 (1|3)(3|4|5|3) 2 ((2|3)(2|4|5|4|5|4|5|2)(3|2|1]+(2|3)((3|5)(5(2|4|1]+6(2|5|1]) 

+(2|5)(3|4|1])(2|1|2|1|2|1|5] - (1|2|1|2)( - 2(3|4)(3|5)(-(2|4|1] 2 + <2|4|1]<2|5|1] + 5(2|5|1] 2 ) 

-(2|5)(3|4)(-(2|4|1] + 2(2|5|1])(3|4|1] - (2|4)(3|5)(-10(2|4|1] - 13<2| 5| 1] )<3| 5| 1] ) [5|4] 

+(1|2)(2(3|4)(3|5)((2|4|1] 2 + 2<2 14| 1] <2| 5 1 1] - 2(2|5|1] 2 ) + (2|5)(3|4)(2(2|4|1] - <2|5|1])<3|4|1] 

+(2|4)(3|5)(5(2|4|1] + 8(2|5|1])(3|5|1]) [5|4|5|4]) }^ 2 (D.4) 



(1|2)(1|5)(2|3)(4|5)[4|3][5|2] 2 2 2 _ (1|2) 2 (1|5) 2 (2|3) 3 [2|1][4|3][5|1][5|2] 
(2|5)(3|4)A {fl ' 2A 3 (2|5) 2 (3|4)(4|5) 

{(4|3|2] 2 (2(4|5) 2 (2|3|2|5) 2 - (1|5)(4|5)(2|3|2|5)(-4(2|5)(4|3|1] + <2|4)(-5<5|2|1] - (5|3|1])) 

+(1|5) 2 (2(2|5) 2 (4|3|1] 2 - (2|4)(2|5)(4|3|1](-5(5|2|1] - (5|3|1]) 

+(2|4) 2 (4(5|2|1] 2 + 3(5|2|1](5|3|1]+(5|3|1] 2 ))) 

+(1|4)(4|3|2]((4|5)(2|3|2|5)((3|5)((4|2|1]+5(4|3|1])+4(3|4)(5|2|1])(5|2|3] 
+(4|5)(5|1|2|3|1|5)(-(2|5)(4|3|1] + (2|4)(-3(5|2|1] - 2(5|3|1])) + (1|5)(5|2|3] 
((2|5)(3|4)(4|3|1](3(5|2|1] + 5(5|3|1]) + (2|4)(3|5)(4|2|1](3(5|2|1] + <5|3|1]) 
+ (2|4)(3|4)(5(5|2|1] 2 + 12(5|2|1](5|3|1] + 3(5|3|1] 2 ))) 

+ (1|4) 2 (((3|5) 2 ((4|2|1] 2 + 3(4|2|1](4|3|1]+4(4|3|1] 2 ) 

-<3|4)<3|5)(-<4|2|1] - 5(4|3|1])(5|2|1] + 2(3|4) 2 (5|2|1] 2 ) (5|2|3] 2 + (4|5) 2 (5|1|3|2|1] 2 
-(4|5)(5|2|1]((3|5)(-2(4|2|1] - 3(4|3|1]) - <3|4)<5|2|1])<5|1|3|2|3]) } M 2 (D.5) 
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